i=1
for(i in 1:length(dependentlist)){
dep=dependentlist[i]
print(dependentlist_eng[i])
load(paste("../data/modeldata/RForigmodeldatawithgeoandgeom_",dep,".RData",sep=""))
preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/RF_geoandgeom/RFwithgeoandgeom_fw_5fold_6p_",dep,"_allpreds",sep=""),kk=1:5,endround = 6,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
print(preds)
importance_ranfor_pset(modeldata=origmodeldata,dependent=dep,pset=5,altdata=origmodeldata,withalt=F)
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=5,kappasum = F,tausum = F,altdata=origmodeldata,withalt=F)
predictors = vector()
for (n in 1:ncol(preds)){
for (n2 in 1:nrow(preds)){
predictors <- c(predictors,as.character(preds[n2,n]))
}
}
print("###################### WITH PREDICTORS from the FW SELECTION ###################")
uniquepredictors <- unique(predictors)
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=uniquepredictors,pset=5,kappasum = F,tausum = F,altdata=origmodeldata,withalt=F)
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"

## [1] "Prediction error at end is: 0.683982683982684"
## [2] "Prediction error at end is: 0.60995670995671"
## [3] "Prediction error at end is: 0.637662337662338"
## [4] "Prediction error at end is: 0.628571428571429"
## [5] "Prediction error at end is: 0.61948051948052"
## [6] "Prediction error at end is: 0.60995670995671"
## k 1
## 1 geom_10m_fl4_L80
## 2 profc_ws7_hr_hr
## 3 geom_10m_fl8_L8
## 4 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr
## 5 vectorruggedness_hr_ws5
## 6 geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr
## k 2
## 1 SGU
## 2 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr
## 3 slope_ws7_10m
## 4 geom_10m_fl3_L22
## 5 geom_10m_fl1_L20
## 6 geom_10m_fl2_L20
## k 3
## 1 geom_10m_fl3_L80
## 2 geom_10m_fl3_L3
## 3 sagaTopographic_Wetness_Index_50m
## 4 dtm_hr_TEXTURE_r10
## 5 crosc_ws3_hr_hr
## 6 geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr
## k 4
## 1 dtm_hr_TEXTURE_r10
## 2 geom_10m_fl3_L50
## 3 geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr
## 4 geom_10m_fl1_L19
## 5 TRI_hr_ws52
## 6 geom_10m_fl1_L25
## k 5
## 1 geom_10m_fl3_L7
## 2 geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr
## 4 Profile_Curvature_50m
## 5 geom_hr_L3_fl1_rplipshape_UE_hr_60cells_hr
## 6 maxic_ws11_10m
## [1] "OBB error with all predictors of allpreds is 0.592592592592593"
## MDA MDG
## geom_10m_fl4_L9 0.002183141 0.4393853
## TRI_hr_ws35 0.002160220 0.3132864
## TRI_hr_ws25 0.001576663 0.2657080
## TRI_hr_ws19 0.001546027 0.2855632
## vectorruggedness_hr_ws13 0.001545295 0.3410383
## vectorruggedness_hr_ws53 0.001417097 0.2337131
## TRI_hr_ws38 0.001398693 0.2275748
## TRI_hr_ws36 0.001301068 0.2962924
## TRI_hr_ws27 0.001292031 0.2284488
## maxic_ws11_hr_hr 0.001283298 0.2312490
## [1] "OBB error with all predictors of allpreds is 0.648148148148148"
## [1] "OOB-error: 0.518518518518518 for predictors geom_10m_fl4_L80 AND profc_ws7_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 1 1 1 1 1.0000000
## 2 0 7 0 4 2 0.4615385
## 3 1 0 2 15 3 0.9047619
## 4 1 3 2 27 5 0.2894737
## 5 0 3 1 12 16 0.5000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.62962962962963"
## [1] "OOB-error: 0.509259259259259 for predictors geom_10m_fl4_L80 AND profc_ws7_hr_hr AND geom_10m_fl8_L8 AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND vectorruggedness_hr_ws5 AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr AND SGU AND slope_ws7_10m AND geom_10m_fl3_L22 AND geom_10m_fl1_L20 AND geom_10m_fl2_L20 AND geom_10m_fl3_L80 AND geom_10m_fl3_L3 AND sagaTopographic_Wetness_Index_50m AND dtm_hr_TEXTURE_r10 AND crosc_ws3_hr_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr AND geom_10m_fl3_L50 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr AND geom_10m_fl1_L19 AND TRI_hr_ws52 AND geom_10m_fl1_L25 AND geom_10m_fl3_L7 AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr AND Profile_Curvature_50m AND geom_hr_L3_fl1_rplipshape_UE_hr_60cells_hr AND maxic_ws11_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 0 3 1 0 1.0000000
## 2 0 4 3 5 1 0.6923077
## 3 0 5 6 8 2 0.7142857
## 4 0 2 3 27 6 0.2894737
## 5 0 0 1 15 16 0.5000000
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"

## [1] "Prediction error at end is: 0.711255411255411"
## [2] "Prediction error at end is: 0.722077922077922"
## [3] "Prediction error at end is: 0.731168831168831"
## [4] "Prediction error at end is: 0.686147186147186"
## [5] "Prediction error at end is: 0.676190476190476"
## [6] "Prediction error at end is: 0.695238095238095"
## k 1
## 1 geom_dtm_10m_hyd_fl5_L6
## 2 vectorruggedness_hr_ws19
## 3 vectorruggedness_hr_ws7
## 4 terraintexture_hr_ws33_tp5
## 5 geom_dtm_10m_hyd_fl5_L37
## 6 geom_10m_fl8_L6
## k 2
## 1 geom_dtm_10m_hyd_fl5_L100
## 2 geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr
## 3 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 4 minic_ws5_hr_hr
## 5 geom_hr_L50m_fl10_rplipedgedensity_UE_hr_5cells_hr
## 6 geom_10m_fl10_L23
## k 3
## 1 geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr
## 2 terraintexture_hr_ws33_tp25
## 3 geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr
## 4 geom_hr_L50m_fl10_rplipshape_UE_hr_5cells_hr
## 5 geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr
## 6 geom_hr_L50m_fl1_rpliprichness_UE_hr_5cells_hr
## k 4
## 1 geom_10m_fl8_L15
## 2 CrossSectionalCurvature_hr
## 3 Relative_Slope_Position_50m
## 4 geom_10m_fl2_L90
## 5 profc_DTM_50m_avg_ws11_50m
## 6 geom_10m_fl2_L100
## k 5
## 1 dtm_hr_TEXTURE_r10
## 2 terraintexture_hr_ws45_tp25
## 3 geom_hr_L50m_fl1_rplippatchdensity_UE_hr_5cells_hr
## 4 geom_10m_fl4_L14
## 5 geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr
## 6 geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr
## [1] "OBB error with all predictors of allpreds is 0.555555555555556"
## MDA MDG
## vectorruggedness_hr_ws7 0.002265043 0.3168125
## vectorruggedness_hr_ws57 0.002261650 0.4070706
## Slope_50m 0.001612188 0.2106449
## vectorruggedness_hr_ws51 0.001359943 0.1927823
## vectorruggedness_hr_ws29 0.001216835 0.1895715
## vectorruggedness_hr_ws59 0.001204661 0.2767824
## terraintexture_hr_ws57_tp25 0.001180399 0.1473424
## vectorruggedness_hr_ws45 0.001155275 0.1984955
## TRI_hr_ws16 0.001153914 0.1778500
## geom_dtm_10m_hyd_fl5_L40.1 0.001148507 0.2754574
## [1] "OBB error with all predictors of allpreds is 0.583333333333333"
## [1] "OOB-error: 0.472222222222222 for predictors geom_dtm_10m_hyd_fl5_L6 AND vectorruggedness_hr_ws19"
## [1] "confusion OOB"
## 2 3 4 5 class.error
## 2 24 3 6 0 0.2727273
## 3 9 10 3 1 0.5652174
## 4 9 6 16 4 0.5428571
## 5 1 3 6 7 0.5882353
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.564814814814815"
## [1] "OOB-error: 0.5 for predictors geom_dtm_10m_hyd_fl5_L6 AND vectorruggedness_hr_ws19 AND vectorruggedness_hr_ws7 AND terraintexture_hr_ws33_tp5 AND geom_dtm_10m_hyd_fl5_L37 AND geom_10m_fl8_L6 AND geom_dtm_10m_hyd_fl5_L100 AND geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND minic_ws5_hr_hr AND geom_hr_L50m_fl10_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl10_L23 AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND terraintexture_hr_ws33_tp25 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr AND geom_hr_L50m_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_5cells_hr AND geom_10m_fl8_L15 AND CrossSectionalCurvature_hr AND Relative_Slope_Position_50m AND geom_10m_fl2_L90 AND profc_DTM_50m_avg_ws11_50m AND geom_10m_fl2_L100 AND dtm_hr_TEXTURE_r10 AND terraintexture_hr_ws45_tp25 AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_5cells_hr AND geom_10m_fl4_L14 AND geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr"
## [1] "confusion OOB"
## 2 3 4 5 class.error
## 2 20 6 7 0 0.3939394
## 3 6 10 7 0 0.5652174
## 4 12 2 19 2 0.4571429
## 5 2 0 10 5 0.7058824
## [1] "1a.3 - Habitat for soil organisms"

## [1] "Prediction error at end is: 0.462337662337662"
## [2] "Prediction error at end is: 0.407359307359307"
## [3] "Prediction error at end is: 0.397835497835498"
## [4] "Prediction error at end is: 0.435497835497835"
## [5] "Prediction error at end is: 0.426406926406926"
## [6] "Prediction error at end is: 0.37012987012987"
## k 1
## 1 TRI_hr_ws14
## 2 NormalizedHeight_10m
## 3 slope_ws11_hr_hr
## 4 TRI_hr_ws5
## 5 geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr
## 6 vectorstrength_hr_ws11_hr
## k 2
## 1 SGU
## 2 profc_DTM_50m_avg_ws3_50m
## 3 geom_10m_fl2_L5
## 4 dtm_hr_TEXTURE_r10
## 5 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr
## k 3
## 1 SGU
## 2 geom_10m_fl3_L33
## 3 slope_DTM_50m_avg_ws5_50m
## 4 Channel_Network_Base_Level_hr
## 5 geom_dtm_10m_hyd_fl5_L5
## 6 geom_10m_fl3_L6
## k 4
## 1 SGU
## 2 minic_DTM_50m_avg_ws9_50m
## 3 geom_10m_fl4_L40
## 4 geom_10m_fl3_L33
## 5 vectorruggedness_hr_ws23
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr
## k 5
## 1 SGU
## 2 elev_ws15_10m
## 3 geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr
## 4 geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr
## 5 geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr
## 6 terraintexture_hr_ws29_tp5
## [1] "OBB error with all predictors of allpreds is 0.388888888888889"
## MDA MDG
## elev_ws11_10m 0.006619780 0.9449497
## elev_ws5_10m 0.006021540 0.7331818
## elev_ws15_10m 0.005803672 0.7382799
## elev_ws3_10m 0.005612781 0.7958222
## Channel_Network_Base_Level_hr 0.004375489 0.6937442
## vectorruggedness_hr_ws5 0.004007954 0.5553209
## vectorruggedness_hr_ws7 0.003688236 0.3846659
## SGU.1 0.003545512 0.7853401
## Normalized_Height_50m 0.003341361 0.5715868
## vectorruggedness_hr_ws9 0.003007679 0.5204331
## [1] "OBB error with all predictors of allpreds is 0.388888888888889"
## [1] "OOB-error: 0.425925925925926 for predictors TRI_hr_ws14 AND NormalizedHeight_10m"
## [1] "confusion OOB"
## 1 2 3 4 class.error
## 1 18 5 1 0 0.2500000
## 2 5 22 2 5 0.3529412
## 3 1 2 15 8 0.4230769
## 4 2 6 9 7 0.7083333
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.37962962962963"
## [1] "OOB-error: 0.314814814814815 for predictors TRI_hr_ws14 AND NormalizedHeight_10m AND slope_ws11_hr_hr AND TRI_hr_ws5 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr AND vectorstrength_hr_ws11_hr AND SGU AND profc_DTM_50m_avg_ws3_50m AND geom_10m_fl2_L5 AND dtm_hr_TEXTURE_r10 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr AND geom_10m_fl3_L33 AND slope_DTM_50m_avg_ws5_50m AND Channel_Network_Base_Level_hr AND geom_dtm_10m_hyd_fl5_L5 AND geom_10m_fl3_L6 AND minic_DTM_50m_avg_ws9_50m AND geom_10m_fl4_L40 AND vectorruggedness_hr_ws23 AND elev_ws15_10m AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr AND terraintexture_hr_ws29_tp5"
## [1] "confusion OOB"
## 1 2 3 4 class.error
## 1 18 5 0 1 0.2500000
## 2 7 25 2 0 0.2647059
## 3 1 2 17 6 0.3461538
## 4 2 1 7 14 0.4166667
## [1] "1a.4 - Habitat for crops"

## [1] "Prediction error at end is: 0.177056277056277"
## [2] "Prediction error at end is: 0.195670995670996"
## [3] "Prediction error at end is: 0.222943722943723"
## [4] "Prediction error at end is: 0.240692640692641"
## [5] "Prediction error at end is: 0.232034632034632"
## [6] "Prediction error at end is: 0.232034632034632"
## k 1
## 1 geom_10m_fl10_L5
## 2 slope_ws3_hr_hr
## 3 geom_10m_fl4_L140
## 4 geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## 6 geom_hr_L3_fl10_rplipedgedensity_UE_hr_5cells_hr
## k 2 k 3
## 1 geom_10m_fl10_L6 geom_10m_fl10_L6
## 2 fischerk_ws27_hr slope_ws5_hr_hr
## 3 fischerk_ws29_hr crosc_DTM_50m_avg_ws5_50m
## 4 planc_DTM_50m_avg_ws3_50m elev_ws11_10m
## 5 PlanCurvature_10m crosc_ws3_hr_hr
## 6 geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr CrossSectionalCurvature_50m
## k 4 k 5
## 1 geom_10m_fl10_L6 slope_ws29_hr_hr
## 2 slope_ws3_hr_hr slope_ws5_hr_hr
## 3 geom_10m_fl8_L23 DiurnalAnisotropicHeating_10m
## 4 geom_10m_fl10_L110 Mid_Slope_Positon_50m
## 5 geom_10m_fl10_L140 minic_ws29_hr_hr
## 6 geom_10m_fl3_L32 TRI_hr_ws8
## [1] "OBB error with all predictors of allpreds is 0.175925925925926"
## MDA MDG
## TRI_hr_ws2 0.010668199 1.1297040
## slope_ws3_hr_hr 0.008896323 0.8565245
## slope_ws5_hr_hr 0.008771970 1.1751347
## TRI_hr_ws1 0.007993862 1.3438221
## MRRTF_hr 0.006752936 0.5735556
## MRVBF_hr 0.005945009 0.7610050
## TRI_hr_ws5 0.005661181 0.6730866
## Slope_hr 0.005382402 0.4765826
## TRI_hr_ws4 0.005264560 0.6865415
## slope_ws19_hr_hr 0.005147809 0.5844093
## [1] "OBB error with all predictors of allpreds is 0.166666666666667"
## [1] "OOB-error: 0.12037037037037 for predictors geom_10m_fl10_L5 AND slope_ws3_hr_hr"
## [1] "confusion OOB"
## 3 4 5 class.error
## 3 4 3 1 0.5000000
## 4 1 59 2 0.0483871
## 5 0 6 32 0.1578947
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.166666666666667"
## [1] "OOB-error: 0.12962962962963 for predictors geom_10m_fl10_L5 AND slope_ws3_hr_hr AND geom_10m_fl4_L140 AND geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl10_L6 AND fischerk_ws27_hr AND fischerk_ws29_hr AND planc_DTM_50m_avg_ws3_50m AND PlanCurvature_10m AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND slope_ws5_hr_hr AND crosc_DTM_50m_avg_ws5_50m AND elev_ws11_10m AND crosc_ws3_hr_hr AND CrossSectionalCurvature_50m AND geom_10m_fl8_L23 AND geom_10m_fl10_L110 AND geom_10m_fl10_L140 AND geom_10m_fl3_L32 AND slope_ws29_hr_hr AND DiurnalAnisotropicHeating_10m AND Mid_Slope_Positon_50m AND minic_ws29_hr_hr AND TRI_hr_ws8"
## [1] "confusion OOB"
## 3 4 5 class.error
## 3 2 6 0 0.7500000
## 4 0 59 3 0.0483871
## 5 0 5 33 0.1315789
## [1] "1c.1 - Average precipitation retention capacity"

## [1] "Prediction error at end is: 0.667099567099567"
## [2] "Prediction error at end is: 0.665367965367965"
## [3] "Prediction error at end is: 0.693506493506494"
## [4] "Prediction error at end is: 0.759307359307359"
## [5] "Prediction error at end is: 0.741125541125541"
## [6] "Prediction error at end is: 0.694372294372294"
## k 1
## 1 geom_10m_fl4_L3
## 2 TRI_hr_ws35
## 3 geom_10m_fl10_L8
## 4 geom_10m_fl10_L14
## 5 geom_10m_fl10_L12
## 6 VerticalDistancetoChannelNetwork_10m
## k 2 k 3
## 1 vectorruggedness_hr_ws55 geom_dtm_10m_hyd_fl5_L6
## 2 Vertical_Distance_to_Channel_Network_50m Channel_Network_Base_Level_hr
## 3 vectorruggedness_hr_ws27 dtm_hr_CONVEX_r5
## 4 elev_ws15_10m geom_dtm_10m_hyd_fl5_L3
## 5 geom_10m_fl10_L6 Convexity_10m
## 6 DiurnalAnisotropicHeating_10m crosc_DTM_50m_avg_ws7_50m
## k 4
## 1 crosc_ws23_hr_hr
## 2 geom_10m_fl2_L4
## 3 terraintexture_hr_ws53_tp5
## 4 slope_ws5_hr_hr
## 5 Plan_Curvature_hr
## 6 minic_ws11_10m
## k 5
## 1 terraintexture_hr_ws37_tp5
## 2 terraintexture_hr_ws41_tp5
## 3 DiurnalAnisotropicHeating_10m
## 4 terraintexture_hr_ws29_tp25
## 5 Relative_Slope_Position_hr
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## [1] "OBB error with all predictors of allpreds is 0.611111111111111"
## MDA MDG
## geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr 0.0020436975 0.40636730
## vectorruggedness_hr_ws57 0.0012897399 0.21516130
## geom_10m_fl3_L110 0.0010810319 0.15107061
## TRI_hr_ws41 0.0010371660 0.14826138
## TRI_hr_ws49 0.0010075752 0.18160740
## vectorruggedness_hr_ws51 0.0009525463 0.18623331
## TRI_hr_ws34 0.0008366258 0.09202898
## vectorruggedness_hr_ws59 0.0008149473 0.27876275
## vectorruggedness_hr_ws41 0.0008117937 0.11315329
## minic_ws7_hr_hr 0.0007986143 0.14062108
## [1] "OBB error with all predictors of allpreds is 0.62962962962963"
## [1] "OOB-error: 0.546296296296296 for predictors geom_10m_fl4_L3 AND TRI_hr_ws35"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 5 14 1 5 1 0.8076923
## 2 4 40 1 0 0 0.1111111
## 3 4 8 0 0 0 1.0000000
## 4 5 9 0 4 0 0.7777778
## 5 3 3 0 1 0 1.0000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.62037037037037"
## [1] "OOB-error: 0.62962962962963 for predictors geom_10m_fl4_L3 AND TRI_hr_ws35 AND geom_10m_fl10_L8 AND geom_10m_fl10_L14 AND geom_10m_fl10_L12 AND VerticalDistancetoChannelNetwork_10m AND vectorruggedness_hr_ws55 AND Vertical_Distance_to_Channel_Network_50m AND vectorruggedness_hr_ws27 AND elev_ws15_10m AND geom_10m_fl10_L6 AND DiurnalAnisotropicHeating_10m AND geom_dtm_10m_hyd_fl5_L6 AND Channel_Network_Base_Level_hr AND dtm_hr_CONVEX_r5 AND geom_dtm_10m_hyd_fl5_L3 AND Convexity_10m AND crosc_DTM_50m_avg_ws7_50m AND crosc_ws23_hr_hr AND geom_10m_fl2_L4 AND terraintexture_hr_ws53_tp5 AND slope_ws5_hr_hr AND Plan_Curvature_hr AND minic_ws11_10m AND terraintexture_hr_ws37_tp5 AND terraintexture_hr_ws41_tp5 AND terraintexture_hr_ws29_tp25 AND Relative_Slope_Position_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 6 16 2 2 0 0.7692308
## 2 6 33 1 5 0 0.2666667
## 3 4 8 0 0 0 1.0000000
## 4 5 12 0 1 0 0.9444444
## 5 4 2 0 1 0 1.0000000
## [1] "1c.1 - Minimum precipitation retention capacity"

## [1] "Prediction error at end is: 0.738528138528138"
## [2] "Prediction error at end is: 0.74025974025974"
## [3] "Prediction error at end is: 0.684848484848485"
## [4] "Prediction error at end is: 0.712554112554112"
## [5] "Prediction error at end is: 0.769264069264069"
## [6] "Prediction error at end is: 0.723376623376623"
## k 1
## 1 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## 2 maxic_ws13_hr_hr
## 3 geom_hr_L50m_fl1_rplippatchnum_UE_hr_60cells_hr
## 4 geom_10m_fl10_L17
## 5 geom_10m_fl1_L5
## 6 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## k 2
## 1 terraintexture_hr_ws53_tp25
## 2 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 3 Vertical_Distance_to_Channel_Network_50m
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr
## 6 geom_hr_L3_fl1_rplippatchdensity_UE_hr_40cells_hr
## k 3
## 1 terraintexture_hr_ws57_tp25
## 2 terraintexture_hr_ws53_tp25
## 3 geom_dtm_10m_hyd_fl5_L30
## 4 geom_hr_L3_fl10_rplipsimpson_UE_hr_5cells_hr
## 5 geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr
## 6 terraintexture_hr_ws17_tp25
## k 4
## 1 profc_ws7_10m
## 2 vectorstrength_hr_ws17_hr
## 3 terraintexture_hr_ws53_tp25
## 4 geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr
## 5 minic_ws15_10m
## 6 terraintexture_hr_ws33_tp25
## k 5
## 1 terraintexture_hr_ws57_tp5
## 2 geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr
## 3 geom_10m_fl10_L9
## 4 terraintexture_hr_ws41_tp5
## 5 geom_10m_fl4_L17
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr
## [1] "OBB error with all predictors of allpreds is 0.694444444444444"
## MDA MDG
## vectorruggedness_hr_ws13 0.0023444153 0.4448103
## TotalCurvature_10m 0.0022660545 0.5706281
## vectorruggedness_hr_ws7 0.0018531709 0.4138484
## vectorruggedness_hr_ws9 0.0016643960 0.3776085
## Channel_Network_Base_Level_hr 0.0012655149 0.3578960
## vectorruggedness_hr_ws15 0.0011921535 0.4308232
## elev_ws5_10m 0.0011401236 0.3988219
## minic_ws11_hr_hr 0.0010589814 0.1553845
## planc_ws29_hr_hr 0.0008921765 0.1827726
## vectorruggedness_hr_ws11 0.0008484432 0.2595412
## [1] "OBB error with all predictors of allpreds is 0.722222222222222"
## [1] "OOB-error: 0.583333333333333 for predictors geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND maxic_ws13_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 4 1 3 3 1 0.6666667
## 2 1 11 8 3 3 0.5769231
## 3 3 9 6 0 4 0.7272727
## 4 3 4 3 15 1 0.4230769
## 5 0 5 2 6 9 0.5909091
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.740740740740741"
## [1] "OOB-error: 0.583333333333333 for predictors geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND maxic_ws13_hr_hr AND geom_hr_L50m_fl1_rplippatchnum_UE_hr_60cells_hr AND geom_10m_fl10_L17 AND geom_10m_fl1_L5 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND terraintexture_hr_ws53_tp25 AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_40cells_hr AND terraintexture_hr_ws57_tp25 AND geom_dtm_10m_hyd_fl5_L30 AND geom_hr_L3_fl10_rplipsimpson_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr AND terraintexture_hr_ws17_tp25 AND profc_ws7_10m AND vectorstrength_hr_ws17_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr AND minic_ws15_10m AND terraintexture_hr_ws33_tp25 AND terraintexture_hr_ws57_tp5 AND geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr AND geom_10m_fl10_L9 AND terraintexture_hr_ws41_tp5 AND geom_10m_fl4_L17"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 2 3 4 3 1.0000000
## 2 0 16 5 4 1 0.3846154
## 3 0 6 10 1 5 0.5454545
## 4 2 5 5 8 6 0.6923077
## 5 1 1 1 8 11 0.5000000
## [1] "1c.2 - Retention capacity for heavy precipitation events"

## [1] "Prediction error at end is: 0.333333333333333"
## [2] "Prediction error at end is: 0.314718614718615"
## [3] "Prediction error at end is: 0.296536796536797"
## [4] "Prediction error at end is: 0.315584415584416"
## [5] "Prediction error at end is: 0.287445887445887"
## [6] "Prediction error at end is: 0.296969696969697"
## k 1
## 1 geom_10m_fl8_L16
## 2 Modified_Catchment_Area_50m
## 3 Catchment_slope_50m
## 4 geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr
## 5 Catchment_Area2_50m
## 6 Total_Curvature_50m
## k 2
## 1 geom_10m_fl1_L8
## 2 crosc_DTM_50m_avg_ws5_50m
## 3 maxic_ws7_hr_hr
## 4 longc_ws3_hr_hr
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_60cells_hr
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr
## k 3
## 1 geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr
## 2 geom_10m_fl2_L8
## 3 crosc_ws11_10m
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr
## 5 Catchment_slope_50m
## 6 SlopeHeight_10m
## k 4 k 5
## 1 geom_10m_fl10_L16 geom_10m_fl3_L11
## 2 geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr geom_dtm_10m_hyd_fl5_L130
## 3 Vertical_Distance_to_Channel_Network_50m geom_10m_fl3_L120
## 4 fischerk_ws53_hr geom_10m_fl1_L9
## 5 FlowLineCurvature_10m Texture_50m
## 6 vectorruggedness_hr_ws35 geom_dtm_10m_hyd_fl5_L110
## [1] "OBB error with all predictors of allpreds is 0.296296296296296"
## MDA MDG
## vectorruggedness_hr_ws31 0.001964531 0.1888182
## crosc_ws13_hr_hr 0.001890218 0.1261239
## vectorruggedness_hr_ws23 0.001786330 0.2351535
## Modified_Catchment_Area_50m 0.001778583 0.2636151
## minic_ws23_hr_hr 0.001735220 0.3290468
## FlowLineCurvature_10m 0.001540143 0.2717212
## crosc_ws15_hr_hr 0.001430740 0.1391839
## vectorruggedness_hr_ws27 0.001424095 0.1962401
## vectorruggedness_hr_ws29 0.001352358 0.1146848
## sagaTopographic_Wetness_Index_50m 0.001347242 0.1522298
## [1] "OBB error with all predictors of allpreds is 0.287037037037037"
## [1] "OOB-error: 0.287037037037037 for predictors geom_10m_fl8_L16 AND Modified_Catchment_Area_50m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 72 2 0 0 1 0.0400000
## 2 7 3 0 0 1 0.7272727
## 3 6 0 0 0 0 1.0000000
## 4 3 2 0 2 1 0.7500000
## 5 6 1 0 1 0 1.0000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.296296296296296"
## [1] "OOB-error: 0.259259259259259 for predictors geom_10m_fl8_L16 AND Modified_Catchment_Area_50m AND Catchment_slope_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND Catchment_Area2_50m AND Total_Curvature_50m AND geom_10m_fl1_L8 AND crosc_DTM_50m_avg_ws5_50m AND maxic_ws7_hr_hr AND longc_ws3_hr_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr AND geom_10m_fl2_L8 AND crosc_ws11_10m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND SlopeHeight_10m AND geom_10m_fl10_L16 AND geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr AND Vertical_Distance_to_Channel_Network_50m AND fischerk_ws53_hr AND FlowLineCurvature_10m AND vectorruggedness_hr_ws35 AND geom_10m_fl3_L11 AND geom_dtm_10m_hyd_fl5_L130 AND geom_10m_fl3_L120 AND geom_10m_fl1_L9 AND Texture_50m AND geom_dtm_10m_hyd_fl5_L110"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 73 1 0 1 0 0.02666667
## 2 8 2 0 1 0 0.81818182
## 3 6 0 0 0 0 1.00000000
## 4 3 0 0 5 0 0.37500000
## 5 7 0 0 1 0 1.00000000
## [1] "1c.3 - groundwater reformation rate"

## [1] "Prediction error at end is: 0.704761904761905"
## [2] "Prediction error at end is: 0.731601731601732"
## [3] "Prediction error at end is: 0.665367965367965"
## [4] "Prediction error at end is: 0.654978354978355"
## [5] "Prediction error at end is: 0.598701298701299"
## [6] "Prediction error at end is: 0.637229437229437"
## k 1
## 1 geom_10m_fl8_L5
## 2 geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl1_rpliprichness_UE_hr_10cells_hr
## 4 geom_10m_fl8_L4
## 5 geom_10m_fl8_L3
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr
## k 2
## 1 SGU
## 2 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 3 slope_DTM_50m_avg_ws7_50m
## 4 SGUT_wTGnew
## 5 maxic_DTM_50m_avg_ws5_50m
## 6 geom_10m_fl10_L6
## k 3
## 1 SGU
## 2 geom_10m_fl8_L50
## 3 CrossSectionalCurvature_hr
## 4 Vertical_Distance_to_Channel_Network_50m
## 5 terraintexture_hr_ws29_tp5
## 6 TotalCurvature_10m
## k 4
## 1 terraintexture_hr_ws49_tp25
## 2 TRI_hr_ws17
## 3 geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr
## 4 fischerk_ws33_hr
## 5 profc_ws7_10m
## 6 terraintexture_hr_ws49_tp5
## k 5
## 1 terraintexture_hr_ws37_tp25
## 2 geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr
## 3 geom_10m_fl8_L7
## 4 vectorstrength_hr_ws13_hr
## 5 geom_10m_fl2_L16
## 6 geom_10m_fl1_L110
## [1] "OBB error with all predictors of allpreds is 0.611111111111111"
## MDA MDG
## elev_ws5_10m 0.001714018 0.3747836
## vectorruggedness_hr_ws29 0.001642631 0.4308794
## vectorruggedness_hr_ws55 0.001265736 0.3826424
## vectorruggedness_hr_ws31 0.001257298 0.3257567
## geom_10m_fl3_L49 0.001248174 0.1705204
## terraintexture_hr_ws53_tp25 0.001155955 0.3237179
## vectorruggedness_hr_ws37 0.001101730 0.2363013
## vectorruggedness_hr_ws59 0.001099157 0.3692933
## vectorruggedness_hr_ws53 0.001097192 0.3506085
## SGU 0.001048640 0.2907233
## [1] "OBB error with all predictors of allpreds is 0.62037037037037"
## [1] "OOB-error: 0.518518518518518 for predictors geom_10m_fl8_L5 AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 4 6 1 0 3 0.7142857
## 2 1 27 3 1 7 0.3076923
## 3 4 5 1 2 2 0.9285714
## 4 0 3 1 7 8 0.6315789
## 5 0 6 3 0 13 0.4090909
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.611111111111111"
## [1] "OOB-error: 0.537037037037037 for predictors geom_10m_fl8_L5 AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl8_L4 AND geom_10m_fl8_L3 AND geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr AND SGU AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND slope_DTM_50m_avg_ws7_50m AND SGUT_wTGnew AND maxic_DTM_50m_avg_ws5_50m AND geom_10m_fl10_L6 AND geom_10m_fl8_L50 AND CrossSectionalCurvature_hr AND Vertical_Distance_to_Channel_Network_50m AND terraintexture_hr_ws29_tp5 AND TotalCurvature_10m AND terraintexture_hr_ws49_tp25 AND TRI_hr_ws17 AND geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr AND fischerk_ws33_hr AND profc_ws7_10m AND terraintexture_hr_ws49_tp5 AND terraintexture_hr_ws37_tp25 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl8_L7 AND vectorstrength_hr_ws13_hr AND geom_10m_fl2_L16 AND geom_10m_fl1_L110"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 3 8 0 0 3 0.7857143
## 2 3 28 3 3 2 0.2820513
## 3 3 4 2 4 1 0.8571429
## 4 0 6 1 8 4 0.5789474
## 5 4 6 0 3 9 0.5909091
## [1] "1c.4 - Potential for providing nutrients for plants"

## [1] "Prediction error at end is: 0.25974025974026"
## [2] "Prediction error at end is: 0.333766233766234"
## [3] "Prediction error at end is: 0.315584415584416"
## [4] "Prediction error at end is: 0.296536796536797"
## [5] "Prediction error at end is: 0.306493506493506"
## [6] "Prediction error at end is: 0.324675324675325"
## k 1
## 1 SGU
## 2 geom_10m_fl10_L3
## 3 vectorstrength_hr_ws17_hr
## 4 geom_10m_fl2_L5
## 5 geom_hr_L50m_fl1_rplipedgedensity_UE_hr_60cells_hr
## 6 VerticalDistancetoChannelNetwork_10m
## k 2
## 1 SGU
## 2 geom_10m_fl8_L11
## 3 geom_hr_L3_fl10_rplippatchdensity_UE_hr_60cells_hr
## 4 geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr
## 5 geom_10m_fl2_L4
## 6 sagaTopographic_Wetness_Index_50m
## k 3
## 1 SGU
## 2 Longitudinal_Curvature_50m
## 3 geom_10m_fl1_L9
## 4 geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr
## 6 geom_10m_fl1_L4
## k 4
## 1 SGU
## 2 ValleyDepth_10m
## 3 VerticalDistancetoChannelNetwork_10m
## 4 Texture_10m
## 5 terraintexture_hr_ws29_tp25
## 6 geom_hr_L50m_fl1_rplippatchdensity_UE_hr_10cells_hr
## k 5
## 1 SGU
## 2 geom_10m_fl3_L12
## 3 Vertical_Distance_to_Channel_Network_hr
## 4 Catchment_slope_50m
## 5 geom_10m_fl1_L30
## 6 geom_10m_fl3_L23
## [1] "OBB error with all predictors of allpreds is 0.305555555555556"
## MDA MDG
## SGU 0.004692489 0.6229274
## SGU.1 0.004663610 0.7811620
## vectorruggedness_hr_ws25 0.004633319 0.6405518
## vectorruggedness_hr_ws29 0.003225298 0.5403808
## vectorruggedness_hr_ws27 0.002638010 0.4342980
## elev_ws3_10m 0.002515565 0.3845359
## elev_ws11_10m 0.002307376 0.4274986
## vectorruggedness_hr_ws7 0.002126472 0.4609184
## minic_ws11_hr_hr 0.002029087 0.2563537
## minic_ws7_hr_hr 0.001951178 0.4521625
## [1] "OBB error with all predictors of allpreds is 0.333333333333333"
## [1] "OOB-error: 0.287037037037037 for predictors SGU AND geom_10m_fl10_L3"
## [1] "confusion OOB"
## 1 3 5 class.error
## 1 55 2 8 0.1538462
## 3 5 0 6 1.0000000
## 5 8 2 22 0.3125000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.287037037037037"
## [1] "OOB-error: 0.222222222222222 for predictors SGU AND geom_10m_fl10_L3 AND vectorstrength_hr_ws17_hr AND geom_10m_fl2_L5 AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_60cells_hr AND VerticalDistancetoChannelNetwork_10m AND geom_10m_fl8_L11 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr AND geom_10m_fl2_L4 AND sagaTopographic_Wetness_Index_50m AND Longitudinal_Curvature_50m AND geom_10m_fl1_L9 AND geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND geom_10m_fl1_L4 AND ValleyDepth_10m AND Texture_10m AND terraintexture_hr_ws29_tp25 AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_10cells_hr AND geom_10m_fl3_L12 AND Vertical_Distance_to_Channel_Network_hr AND Catchment_slope_50m AND geom_10m_fl1_L30 AND geom_10m_fl3_L23"
## [1] "confusion OOB"
## 1 3 5 class.error
## 1 60 1 4 0.07692308
## 3 7 0 4 1.00000000
## 5 8 0 24 0.25000000
## [1] "1c.5 - Potential as a CO2 sink"

## [1] "Prediction error at end is: 0.541991341991342"
## [2] "Prediction error at end is: 0.466666666666667"
## [3] "Prediction error at end is: 0.42987012987013"
## [4] "Prediction error at end is: 0.382251082251082"
## [5] "Prediction error at end is: 0.410822510822511"
## [6] "Prediction error at end is: 0.392640692640693"
## k 1
## 1 slope_ws19_hr_hr
## 2 fischerk_ws19_hr
## 3 vectorstrength_hr_ws27_hr
## 4 geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr
## 5 minic_DTM_50m_avg_ws11_50m
## 6 TRI_hr_ws18
## k 2
## 1 slope_ws19_hr_hr
## 2 maxic_DTM_50m_avg_ws5_50m
## 3 vectorruggedness_hr_ws7
## 4 geom_hr_L3_fl1_rplipshape_UE_hr_40cells_hr
## 5 geom_10m_fl3_L110
## 6 geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr
## k 3
## 1 CrossSectionalCurvature_hr
## 2 TotalCurvature_10m
## 3 Topographic_Wetness_Index_50m
## 4 Channel_Network_Base_Level_hr
## 5 geom_hr_L50m_fl1_rplipdominance_UE_hr_20cells_hr
## 6 longc_ws5_hr_hr
## k 4 k 5
## 1 crosc_ws5_hr_hr SGU
## 2 slope_ws23_hr_hr vectorruggedness_hr_ws35
## 3 geom_dtm_10m_hyd_fl5_L29 MassBalanceIndex_10m
## 4 geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr crosc_ws3_10m
## 5 CrossSectionalCurvature_hr SGUT_wTGnew
## 6 geom_10m_fl2_L140 geom_10m_fl4_L19
## [1] "OBB error with all predictors of allpreds is 0.342592592592593"
## MDA MDG
## elev_ws5_10m 0.006335650 0.8492952
## elev_ws3_10m 0.005769299 0.8009306
## vectorruggedness_hr_ws29 0.005304706 0.6185339
## vectorruggedness_hr_ws9 0.005133086 0.5318235
## vectorruggedness_hr_ws15 0.004717844 0.6522424
## elev_ws11_10m 0.004677776 0.7129504
## elev_ws15_10m 0.004638132 0.7054063
## vectorruggedness_hr_ws7 0.004200919 0.5919718
## Channel_Network_Base_Level_hr 0.004171763 0.5533772
## vectorruggedness_hr_ws41 0.004047411 0.5727091
## [1] "OBB error with all predictors of allpreds is 0.324074074074074"
## [1] "OOB-error: 0.324074074074074 for predictors slope_ws19_hr_hr AND fischerk_ws19_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 48 0 2 1 0 0.05882353
## 2 0 0 1 0 0 1.00000000
## 3 4 1 3 10 0 0.83333333
## 4 4 0 7 22 0 0.33333333
## 5 1 0 1 3 0 1.00000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.324074074074074"
## [1] "OOB-error: 0.259259259259259 for predictors slope_ws19_hr_hr AND fischerk_ws19_hr AND vectorstrength_hr_ws27_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr AND minic_DTM_50m_avg_ws11_50m AND TRI_hr_ws18 AND maxic_DTM_50m_avg_ws5_50m AND vectorruggedness_hr_ws7 AND geom_hr_L3_fl1_rplipshape_UE_hr_40cells_hr AND geom_10m_fl3_L110 AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND CrossSectionalCurvature_hr AND TotalCurvature_10m AND Topographic_Wetness_Index_50m AND Channel_Network_Base_Level_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_20cells_hr AND longc_ws5_hr_hr AND crosc_ws5_hr_hr AND slope_ws23_hr_hr AND geom_dtm_10m_hyd_fl5_L29 AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND geom_10m_fl2_L140 AND SGU AND vectorruggedness_hr_ws35 AND MassBalanceIndex_10m AND crosc_ws3_10m AND SGUT_wTGnew AND geom_10m_fl4_L19"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 45 0 2 4 0 0.1176471
## 2 0 0 1 0 0 1.0000000
## 3 2 0 9 7 0 0.5000000
## 4 5 0 2 26 0 0.2121212
## 5 1 0 0 4 0 1.0000000
## [1] "1d.1 - Potential for retention of heavy metals"

## [1] "Prediction error at end is: 0.444588744588745"
## [2] "Prediction error at end is: 0.481385281385281"
## [3] "Prediction error at end is: 0.407792207792208"
## [4] "Prediction error at end is: 0.398268398268398"
## [5] "Prediction error at end is: 0.361471861471861"
## [6] "Prediction error at end is: 0.38961038961039"
## k 1
## 1 SGU
## 2 vectorruggedness_hr_ws45
## 3 Channel_Network_Base_Level_hr
## 4 Catchment_slope_50m
## 5 geom_10m_fl4_L47
## 6 fischerk_ws17_hr
## k 2
## 1 geom_10m_fl8_L60
## 2 Convexity_50m
## 3 geom_10m_fl8_L50
## 4 Relative_Slope_Position_50m
## 5 geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr
## 6 geom_10m_fl4_L14
## k 3
## 1 geom_10m_fl8_L8
## 2 geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr
## 3 geom_dtm_10m_hyd_fl5_L60
## 4 geom_10m_fl4_L23
## 5 geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr
## 6 Valley_Depth_hr
## k 4
## 1 geom_10m_fl8_L12
## 2 geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr
## 3 geom_10m_fl3_L25
## 4 profc_DTM_50m_avg_ws7_50m
## 5 geom_10m_fl1_L50
## 6 geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr
## k 5
## 1 SGU
## 2 geom_10m_fl1_L42
## 3 TangentialCurvature_10m
## 4 ChannelNetworkBaseLevel_10m
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr
## 6 vectorruggedness_hr_ws59
## [1] "OBB error with all predictors of allpreds is 0.324074074074074"
## MDA MDG
## elev_ws3_10m 0.004953419 0.7750291
## vectorruggedness_hr_ws27 0.004922231 0.7177749
## elev_ws15_10m 0.003974022 0.5810815
## vectorruggedness_hr_ws31 0.003474228 0.5501727
## elev_ws11_10m 0.003352464 0.5704524
## vectorruggedness_hr_ws23 0.003267820 0.4449253
## vectorruggedness_hr_ws41 0.003247652 0.6457645
## elev_ws5_10m 0.003193369 0.4353484
## vectorruggedness_hr_ws43 0.003028918 0.3882980
## vectorruggedness_hr_ws25 0.002957736 0.3383449
## [1] "OBB error with all predictors of allpreds is 0.342592592592593"
## [1] "OOB-error: 0.305555555555556 for predictors SGU AND vectorruggedness_hr_ws45"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 36 1 1 0 4 0.1428571
## 2 2 2 0 1 2 0.7142857
## 3 6 0 2 0 2 0.8000000
## 4 2 1 0 2 2 0.7142857
## 5 7 1 1 0 33 0.2142857
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.324074074074074"
## [1] "OOB-error: 0.277777777777778 for predictors SGU AND vectorruggedness_hr_ws45 AND Channel_Network_Base_Level_hr AND Catchment_slope_50m AND geom_10m_fl4_L47 AND fischerk_ws17_hr AND geom_10m_fl8_L60 AND Convexity_50m AND geom_10m_fl8_L50 AND Relative_Slope_Position_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr AND geom_10m_fl4_L14 AND geom_10m_fl8_L8 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND geom_dtm_10m_hyd_fl5_L60 AND geom_10m_fl4_L23 AND geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr AND Valley_Depth_hr AND geom_10m_fl8_L12 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND geom_10m_fl3_L25 AND profc_DTM_50m_avg_ws7_50m AND geom_10m_fl1_L50 AND geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr AND geom_10m_fl1_L42 AND TangentialCurvature_10m AND ChannelNetworkBaseLevel_10m AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND vectorruggedness_hr_ws59"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 38 0 1 0 3 0.0952381
## 2 3 1 0 0 3 0.8571429
## 3 5 0 2 0 3 0.8000000
## 4 2 0 0 0 5 1.0000000
## 5 5 0 0 0 37 0.1190476
## [1] "1d.2 - Potential for transforming organic contaminants"

## [1] "Prediction error at end is: 0.6"
## [2] "Prediction error at end is: 0.572294372294372"
## [3] "Prediction error at end is: 0.654978354978355"
## [4] "Prediction error at end is: 0.617748917748918"
## [5] "Prediction error at end is: 0.627705627705628"
## [6] "Prediction error at end is: 0.637662337662338"
## k 1
## 1 geom_10m_fl10_L21
## 2 geom_dtm_10m_hyd_fl5_L110
## 3 geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr
## 4 minic_DTM_50m_avg_ws11_50m
## 5 geom_10m_fl10_L150
## 6 maxic_DTM_50m_avg_ws11_50m
## k 2 k 3
## 1 geom_dtm_10m_hyd_fl5_L8 geom_10m_fl10_L90
## 2 geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr elev_ws3_10m
## 3 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr planc_DTM_50m_avg_ws3_50m
## 4 geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr geom_10m_fl10_L120
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr maxic_ws3_hr_hr
## 6 geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr geom_10m_fl8_L120
## k 4
## 1 geom_10m_fl10_L22
## 2 PlanCurvature_10m
## 3 geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr
## 4 geom_10m_fl3_L6
## 5 geom_10m_fl10_L17
## 6 geom_10m_fl10_L19
## k 5
## 1 geom_10m_fl10_L22
## 2 elev_ws15_10m
## 3 geom_10m_fl8_L9
## 4 Convexity_50m
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr
## 6 minic_ws3_10m
## [1] "OBB error with all predictors of allpreds is 0.546296296296296"
## MDA MDG
## elev_ws5_10m 0.003169042 0.5935770
## elev_ws15_10m 0.002780281 0.7086013
## elev_ws11_10m 0.002493002 0.5607704
## vectorruggedness_hr_ws25 0.002304768 0.5029180
## vectorruggedness_hr_ws41 0.002271230 0.3682756
## vectorruggedness_hr_ws31 0.002191556 0.3236921
## vectorruggedness_hr_ws9 0.002137198 0.3201390
## elev_ws3_10m 0.001802644 0.4526344
## vectorruggedness_hr_ws45 0.001653024 0.3398623
## TotalCurvature_10m 0.001575140 0.4132897
## [1] "OBB error with all predictors of allpreds is 0.537037037037037"
## [1] "OOB-error: 0.490740740740741 for predictors geom_10m_fl10_L21 AND geom_dtm_10m_hyd_fl5_L110"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 3 4 2 0 1 0.7000000
## 2 1 14 4 0 3 0.3636364
## 3 1 7 11 0 9 0.6071429
## 4 0 2 0 0 7 1.0000000
## 5 1 5 6 0 27 0.3076923
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.546296296296296"
## [1] "OOB-error: 0.453703703703704 for predictors geom_10m_fl10_L21 AND geom_dtm_10m_hyd_fl5_L110 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr AND minic_DTM_50m_avg_ws11_50m AND geom_10m_fl10_L150 AND maxic_DTM_50m_avg_ws11_50m AND geom_dtm_10m_hyd_fl5_L8 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl10_L90 AND elev_ws3_10m AND planc_DTM_50m_avg_ws3_50m AND geom_10m_fl10_L120 AND maxic_ws3_hr_hr AND geom_10m_fl8_L120 AND geom_10m_fl10_L22 AND PlanCurvature_10m AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr AND geom_10m_fl3_L6 AND geom_10m_fl10_L17 AND geom_10m_fl10_L19 AND elev_ws15_10m AND geom_10m_fl8_L9 AND Convexity_50m AND geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr AND minic_ws3_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 2 3 4 0 1 0.8000000
## 2 3 10 7 0 2 0.5454545
## 3 0 5 17 0 6 0.3928571
## 4 0 0 2 0 7 1.0000000
## 5 0 3 5 1 30 0.2307692
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"

## [1] "Prediction error at end is: 0.0826839826839827"
## [2] "Prediction error at end is: 0.0463203463203463"
## [3] "Prediction error at end is: 0.0463203463203463"
## [4] "Prediction error at end is: 0.0372294372294372"
## [5] "Prediction error at end is: 0.0372294372294372"
## [6] "Prediction error at end is: 0.0372294372294372"
## k 1
## 1 Slope_Height_50m
## 2 geom_10m_fl2_L15
## 3 geom_10m_fl2_L16
## 4 geom_10m_fl2_L17
## 5 geom_10m_fl2_L18
## 6 geom_hr_L3_fl10_rplipedgedensity_UE_hr_40cells_hr
## k 2
## 1 geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr
## 2 planc_DTM_50m_avg_ws5_50m
## 3 geom_10m_fl2_L19
## 4 geom_10m_fl2_L22
## 5 geom_10m_fl4_L3
## 6 geom_10m_fl8_L17
## k 3
## 1 planc_ws19_hr_hr
## 2 MinimalCurvature_10m
## 3 crosc_ws11_hr_hr
## 4 geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr
## 5 geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr
## 6 fischerk_ws31_hr
## k 4
## 1 maxic_ws7_hr_hr
## 2 geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr
## 3 geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr
## 4 geom_10m_fl10_L3
## 5 minic_ws15_hr_hr
## 6 geom_10m_fl10_L11
## k 5
## 1 geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr
## 2 geom_10m_fl2_L19
## 3 geom_10m_fl4_L8
## 4 geom_10m_fl1_L60
## 5 geom_10m_fl2_L60
## 6 minic_ws7_10m
## [1] "OBB error with all predictors of allpreds is 0.037037037037037"
## MDA MDG
## geom_10m_fl4_L46.1 0.0003617464 0.024221704
## Slope_Height_50m 0.0003543389 0.126936116
## crosc_DTM_50m_avg_ws7_50m 0.0003414634 0.006857143
## geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr 0.0002980622 0.045201805
## Slope_hr 0.0002887207 0.015066667
## geom_10m_fl2_L18 0.0002676692 0.041457817
## maxic_ws15_10m 0.0002638353 0.005047978
## geom_10m_fl4_L39.1 0.0002522648 0.036213276
## geom_10m_fl10_L13.1 0.0002439024 0.008974359
## geom_10m_fl3_L16 0.0002380952 0.010909091
## [1] "OBB error with all predictors of allpreds is 0.037037037037037"
## [1] "OOB-error: 0.0185185185185185 for predictors Slope_Height_50m AND geom_10m_fl2_L15"
## [1] "confusion OOB"
## 4 5 class.error
## 4 2 2 0.5
## 5 0 104 0.0
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.037037037037037"
## [1] "OOB-error: 0.0277777777777778 for predictors Slope_Height_50m AND geom_10m_fl2_L15 AND geom_10m_fl2_L16 AND geom_10m_fl2_L17 AND geom_10m_fl2_L18 AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_40cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND planc_DTM_50m_avg_ws5_50m AND geom_10m_fl2_L19 AND geom_10m_fl2_L22 AND geom_10m_fl4_L3 AND geom_10m_fl8_L17 AND planc_ws19_hr_hr AND MinimalCurvature_10m AND crosc_ws11_hr_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr AND fischerk_ws31_hr AND maxic_ws7_hr_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND geom_10m_fl10_L3 AND minic_ws15_hr_hr AND geom_10m_fl10_L11 AND geom_10m_fl4_L8 AND geom_10m_fl1_L60 AND geom_10m_fl2_L60 AND minic_ws7_10m"
## [1] "confusion OOB"
## 4 5 class.error
## 4 1 3 0.75
## 5 0 104 0.00
## [1] "1d.4 - Potential for retention of water-soluble contaminants"

## [1] "Prediction error at end is: 0.583549783549784"
## [2] "Prediction error at end is: 0.582683982683983"
## [3] "Prediction error at end is: 0.601298701298701"
## [4] "Prediction error at end is: 0.583116883116883"
## [5] "Prediction error at end is: 0.582683982683983"
## [6] "Prediction error at end is: 0.564069264069264"
## k 1
## 1 geom_10m_fl8_L5
## 2 Relative_Slope_Position_50m
## 3 minic_DTM_50m_avg_ws5_50m
## 4 geom_10m_fl2_L15
## 5 geom_10m_fl1_L7
## 6 Plan_Curvature_hr
## k 2
## 1 geom_10m_fl3_L120
## 2 geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr
## 4 geom_10m_fl3_L130
## 5 geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr
## 6 geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr
## k 3 k 4
## 1 geom_10m_fl8_L4 geom_10m_fl8_L110
## 2 Relative_Slope_Position_50m longc_DTM_50m_avg_ws11_50m
## 3 TRI_hr_ws53 vectorruggedness_hr_ws35
## 4 maxic_ws5_hr_hr slope_ws15_10m
## 5 geom_10m_fl8_L5 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 6 NormalizedHeight_10m profc_ws5_hr_hr
## k 5
## 1 geom_10m_fl10_L14
## 2 geom_10m_fl10_L10
## 3 geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr
## 5 geom_dtm_10m_hyd_fl5_L32
## 6 geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr
## [1] "OBB error with all predictors of allpreds is 0.518518518518518"
## MDA MDG
## sagaTopographic_Wetness_Index_50m 0.002988331 0.3318349
## vectorruggedness_hr_ws39 0.002352883 0.3790400
## vectorruggedness_hr_ws27 0.002240995 0.5211155
## vectorruggedness_hr_ws33 0.001947569 0.3315818
## crosc_ws29_hr_hr 0.001803880 0.2807113
## vectorruggedness_hr_ws35 0.001673071 0.2451025
## vectorruggedness_hr_ws25 0.001626234 0.3638782
## vectorruggedness_hr_ws23 0.001529754 0.3604870
## vectorruggedness_hr_ws53 0.001490791 0.3899991
## elev_ws15_10m 0.001467265 0.2676537
## [1] "OBB error with all predictors of allpreds is 0.509259259259259"
## [1] "OOB-error: 0.425925925925926 for predictors geom_10m_fl8_L5 AND Relative_Slope_Position_50m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 45 0 0 2 0 0.04255319
## 2 6 0 2 6 0 1.00000000
## 3 5 2 5 5 0 0.70588235
## 4 5 3 2 12 0 0.45454545
## 5 1 0 1 6 0 1.00000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.509259259259259"
## [1] "OOB-error: 0.481481481481481 for predictors geom_10m_fl8_L5 AND Relative_Slope_Position_50m AND minic_DTM_50m_avg_ws5_50m AND geom_10m_fl2_L15 AND geom_10m_fl1_L7 AND Plan_Curvature_hr AND geom_10m_fl3_L120 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr AND geom_10m_fl3_L130 AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl8_L4 AND TRI_hr_ws53 AND maxic_ws5_hr_hr AND NormalizedHeight_10m AND geom_10m_fl8_L110 AND longc_DTM_50m_avg_ws11_50m AND vectorruggedness_hr_ws35 AND slope_ws15_10m AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND profc_ws5_hr_hr AND geom_10m_fl10_L14 AND geom_10m_fl10_L10 AND geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND geom_dtm_10m_hyd_fl5_L32 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 42 2 3 0 0 0.1063830
## 2 9 0 1 4 0 1.0000000
## 3 5 1 4 6 1 0.7647059
## 4 7 2 2 10 1 0.5454545
## 5 0 0 3 5 0 1.0000000
## [1] "1d.5 - Potential as buffer for acidic contaminants"

## [1] "Prediction error at end is: 0.676190476190476"
## [2] "Prediction error at end is: 0.722077922077922"
## [3] "Prediction error at end is: 0.712987012987013"
## [4] "Prediction error at end is: 0.703030303030303"
## [5] "Prediction error at end is: 0.667099567099567"
## [6] "Prediction error at end is: 0.603030303030303"
## k 1 k 2
## 1 geom_dtm_10m_hyd_fl5_L110 geom_dtm_10m_hyd_fl5_L140
## 2 Catchment_area_50m planc_DTM_50m_avg_ws9_50m
## 3 geom_10m_fl4_L40 SGU
## 4 planc_ws5_hr_hr crosc_ws5_hr_hr
## 5 geom_dtm_10m_hyd_fl5_L42 geom_10m_fl10_L3
## 6 terraintexture_hr_ws53_tp25 geom_dtm_10m_hyd_fl5_L4
## k 3
## 1 geom_hr_L3_fl10_rplipshannon_UE_hr_10cells_hr
## 2 geom_dtm_10m_hyd_fl5_L110
## 3 maxic_ws3_hr_hr
## 4 terraintexture_hr_ws41_tp5
## 5 planc_ws29_hr_hr
## 6 LS_Factor_50m
## k 4 k 5
## 1 geom_10m_fl1_L12 geom_10m_fl1_L8
## 2 geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr geom_10m_fl4_L49
## 3 StandardizedHeight_10m elev_ws5_10m
## 4 geom_dtm_10m_hyd_fl5_L5 Catchment_slope_50m
## 5 geom_10m_fl1_L4 slope_ws29_hr_hr
## 6 geom_hr_L50m_fl10_rpliprichness_UE_hr_5cells_hr FlowLineCurvature_10m
## [1] "OBB error with all predictors of allpreds is 0.546296296296296"
## MDA MDG
## Standardized_Height_50m 0.003973304 0.8203215
## Vertical_Distance_to_Channel_Network_50m 0.003512193 0.4564479
## Modified_Catchment_Area_50m 0.003471425 0.6536779
## Normalized_Height_50m 0.003022192 0.6182740
## Catchment_area_50m 0.002489700 0.4133209
## Catchment_Area2_50m 0.002398727 0.4977011
## SGU.1 0.002150208 0.5578115
## geom_10m_fl10_L3 0.001996288 0.2505707
## SGUT_wTGnew.1 0.001974247 0.4986677
## vectorruggedness_hr_ws19 0.001684300 0.2947276
## [1] "OBB error with all predictors of allpreds is 0.592592592592593"
## [1] "OOB-error: 0.574074074074074 for predictors geom_dtm_10m_hyd_fl5_L110 AND Catchment_area_50m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 25 9 4 1 0 0.3589744
## 2 10 8 4 0 2 0.6666667
## 3 6 3 8 2 1 0.6000000
## 4 2 4 2 5 3 0.6875000
## 5 1 2 3 3 0 1.0000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.592592592592593"
## [1] "OOB-error: 0.388888888888889 for predictors geom_dtm_10m_hyd_fl5_L110 AND Catchment_area_50m AND geom_10m_fl4_L40 AND planc_ws5_hr_hr AND geom_dtm_10m_hyd_fl5_L42 AND terraintexture_hr_ws53_tp25 AND geom_dtm_10m_hyd_fl5_L140 AND planc_DTM_50m_avg_ws9_50m AND SGU AND crosc_ws5_hr_hr AND geom_10m_fl10_L3 AND geom_dtm_10m_hyd_fl5_L4 AND geom_hr_L3_fl10_rplipshannon_UE_hr_10cells_hr AND maxic_ws3_hr_hr AND terraintexture_hr_ws41_tp5 AND planc_ws29_hr_hr AND LS_Factor_50m AND geom_10m_fl1_L12 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr AND StandardizedHeight_10m AND geom_dtm_10m_hyd_fl5_L5 AND geom_10m_fl1_L4 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_5cells_hr AND geom_10m_fl1_L8 AND geom_10m_fl4_L49 AND elev_ws5_10m AND Catchment_slope_50m AND slope_ws29_hr_hr AND FlowLineCurvature_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 31 5 3 0 0 0.2051282
## 2 6 16 1 1 0 0.3333333
## 3 8 1 9 2 0 0.5500000
## 4 1 2 4 8 1 0.5000000
## 5 1 1 2 3 2 0.7777778
i=1
for(i in 1:length(dependentlist)){
dep=dependentlist[i]
print(dependentlist_eng[i])
load(paste("../data/modeldata/RForigmodeldatawithgeoandgeom_",dep,".RData",sep=""))
preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/RF_localterrain/RFwithgeoandgeom_fw_5fold_6p_",dep,"_localterrain",sep=""),kk=1:5,endround = 6,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
print(preds)
importance_ranfor_pset(modeldata=origmodeldata,dependent=dep,pset=1,altdata=origmodeldata,withalt=F)
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=1,kappasum = F,tausum = F,altdata=origmodeldata,withalt=F)
predictors = vector()
for (n in 1:ncol(preds)){
for (n2 in 1:nrow(preds)){
predictors <- c(predictors,as.character(preds[n2,n]))
}
}
print("###################### WITH PREDICTORS from the FW SELECTION ###################")
uniquepredictors <- unique(predictors)
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=uniquepredictors,pset=1,kappasum = F,tausum = F,altdata=origmodeldata,withalt=F)
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"

## [1] "Prediction error at end is: 0.703463203463203"
## [2] "Prediction error at end is: 0.731168831168831"
## [3] "Prediction error at end is: 0.704329004329004"
## [4] "Prediction error at end is: 0.686580086580087"
## [5] "Prediction error at end is: 0.667532467532468"
## [6] "Prediction error at end is: 0.685714285714286"
## k 1 k 2
## 1 Slope_50m minic_ws23_hr_hr
## 2 CrossSectionalCurvature_hr profc_ws13_hr_hr
## 3 longc_DTM_50m_avg_ws3_50m slope_DTM_50m_avg_ws11_50m
## 4 MaximalCurvature_10m Convexity_50m
## 5 longc_DTM_50m_avg_ws5_50m minic_ws19_hr_hr
## 6 profc_DTM_50m_avg_ws3_50m profc_ws11_hr_hr
## k 3 k 4
## 1 maxic_ws9_hr_hr crosc_ws3_10m
## 2 crosc_ws7_10m Total_Curvature_50m
## 3 slope_DTM_50m_avg_ws11_50m Plan_Curvature_hr
## 4 planc_ws5_10m planc_DTM_50m_avg_ws3_50m
## 5 maxic_ws5_hr_hr planc_ws3_hr_hr
## 6 minic_ws11_hr_hr Slope_50m
## k 5
## 1 crosc_ws19_hr_hr
## 2 Plan_Curvature_hr
## 3 longc_ws5_hr_hr
## 4 crosc_ws23_hr_hr
## 5 maxic_ws3_hr_hr
## 6 Maximal_Curvature_hr
## [1] "OBB error with all predictors of localterrain is 0.574074074074074"
## MDA MDG
## slope_ws3_10m 0.004583860 1.0882822
## slope_ws29_hr_hr 0.003768368 0.6986251
## Slope_50m 0.003159166 1.0467980
## Slope_10m 0.003133764 0.6413639
## slope_ws11_hr_hr 0.003024752 0.7680144
## slope_ws5_10m 0.002742490 1.1343291
## slope_ws7_10m 0.002537431 0.6957050
## maxic_ws3_10m 0.002441898 0.6120702
## slope_DTM_50m_avg_ws9_50m 0.002395678 0.8567379
## slope_ws23_hr_hr 0.002221673 0.7688141
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## [1] "OOB-error: 0.564814814814815 for predictors Slope_50m AND CrossSectionalCurvature_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 1 1 2 0 1.0000000
## 2 0 3 4 6 0 0.7692308
## 3 0 3 7 7 4 0.6666667
## 4 0 4 8 19 7 0.5000000
## 5 0 1 3 10 18 0.4375000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.648148148148148"
## [1] "OOB-error: 0.611111111111111 for predictors Slope_50m AND CrossSectionalCurvature_hr AND longc_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND longc_DTM_50m_avg_ws5_50m AND profc_DTM_50m_avg_ws3_50m AND minic_ws23_hr_hr AND profc_ws13_hr_hr AND slope_DTM_50m_avg_ws11_50m AND Convexity_50m AND minic_ws19_hr_hr AND profc_ws11_hr_hr AND maxic_ws9_hr_hr AND crosc_ws7_10m AND planc_ws5_10m AND maxic_ws5_hr_hr AND minic_ws11_hr_hr AND crosc_ws3_10m AND Total_Curvature_50m AND Plan_Curvature_hr AND planc_DTM_50m_avg_ws3_50m AND planc_ws3_hr_hr AND crosc_ws19_hr_hr AND longc_ws5_hr_hr AND crosc_ws23_hr_hr AND maxic_ws3_hr_hr AND Maximal_Curvature_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 2 1 1 0 1.0000000
## 2 0 1 1 11 0 0.9230769
## 3 0 1 3 14 3 0.8571429
## 4 0 1 9 20 8 0.4736842
## 5 0 0 2 12 18 0.4375000
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"

## [1] "Prediction error at end is: 0.769264069264069"
## [2] "Prediction error at end is: 0.546320346320346"
## [3] "Prediction error at end is: 0.564502164502165"
## [4] "Prediction error at end is: 0.584415584415584"
## [5] "Prediction error at end is: 0.584415584415584"
## [6] "Prediction error at end is: 0.583982683982684"
## k 1 k 2
## 1 Slope_50m Plan_Curvature_hr
## 2 GeneralCurvature_10m crosc_ws13_hr_hr
## 3 Flow_Line_Curvature_hr maxic_ws13_hr_hr
## 4 LongitudinalCurvature_10m slope_DTM_50m_avg_ws3_50m
## 5 Minimal_Curvature_hr slope_ws29_hr_hr
## 6 longc_ws5_10m Profile_Curvature_50m
## k 3 k 4
## 1 Total_Curvature_50m maxic_ws9_hr_hr
## 2 Convexity_50m slope_DTM_50m_avg_ws5_50m
## 3 maxic_ws11_10m maxic_ws5_hr_hr
## 4 slope_DTM_50m_avg_ws3_50m maxic_ws19_hr_hr
## 5 profc_ws15_10m maxic_ws13_hr_hr
## 6 profc_ws11_10m planc_ws7_10m
## k 5
## 1 crosc_ws13_hr_hr
## 2 maxic_ws29_hr_hr
## 3 maxic_ws9_hr_hr
## 4 planc_ws9_hr_hr
## 5 Tangential_Curvature_hr
## 6 slope_ws29_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.509259259259259"
## MDA MDG
## Slope_50m 0.006717291 1.1720831
## Slope_10m 0.004280599 0.8046896
## slope_ws19_hr_hr 0.003465982 0.6774488
## slope_DTM_50m_avg_ws9_50m 0.003410984 0.8973186
## maxic_ws13_hr_hr 0.003357477 0.6163522
## slope_ws15_10m 0.003333417 0.7061658
## maxic_ws9_hr_hr 0.003324579 0.8014618
## MaximalCurvature_10m 0.003299902 0.9929497
## slope_ws11_10m 0.003259476 0.7366998
## slope_ws3_10m 0.003211550 0.8055028
## [1] "OBB error with all predictors of localterrain is 0.537037037037037"
## [1] "OOB-error: 0.481481481481481 for predictors Slope_50m AND GeneralCurvature_10m"
## [1] "confusion OOB"
## 2 3 4 5 class.error
## 2 22 7 4 0 0.3333333
## 3 6 8 7 2 0.6521739
## 4 5 5 19 6 0.4571429
## 5 1 2 7 7 0.5882353
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.555555555555556"
## [1] "OOB-error: 0.444444444444444 for predictors Slope_50m AND GeneralCurvature_10m AND Flow_Line_Curvature_hr AND LongitudinalCurvature_10m AND Minimal_Curvature_hr AND longc_ws5_10m AND Plan_Curvature_hr AND crosc_ws13_hr_hr AND maxic_ws13_hr_hr AND slope_DTM_50m_avg_ws3_50m AND slope_ws29_hr_hr AND Profile_Curvature_50m AND Total_Curvature_50m AND Convexity_50m AND maxic_ws11_10m AND profc_ws15_10m AND profc_ws11_10m AND maxic_ws9_hr_hr AND slope_DTM_50m_avg_ws5_50m AND maxic_ws5_hr_hr AND maxic_ws19_hr_hr AND planc_ws7_10m AND maxic_ws29_hr_hr AND planc_ws9_hr_hr AND Tangential_Curvature_hr"
## [1] "confusion OOB"
## 2 3 4 5 class.error
## 2 21 5 7 0 0.3636364
## 3 5 10 7 1 0.5652174
## 4 6 4 23 2 0.3428571
## 5 1 1 9 6 0.6470588
## [1] "1a.3 - Habitat for soil organisms"

## [1] "Prediction error at end is: 0.675324675324675"
## [2] "Prediction error at end is: 0.694372294372294"
## [3] "Prediction error at end is: 0.648917748917749"
## [4] "Prediction error at end is: 0.602597402597403"
## [5] "Prediction error at end is: 0.528571428571429"
## [6] "Prediction error at end is: 0.567099567099567"
## k 1 k 2
## 1 minic_ws9_hr_hr minic_ws9_hr_hr
## 2 slope_ws19_hr_hr profc_DTM_50m_avg_ws3_50m
## 3 minic_DTM_50m_avg_ws7_50m crosc_ws13_hr_hr
## 4 slope_DTM_50m_avg_ws9_50m dtm_hr_CONVEX_r30
## 5 minic_ws5_hr_hr planc_ws15_hr_hr
## 6 planc_ws9_hr_hr crosc_DTM_50m_avg_ws11_50m
## k 3 k 4
## 1 crosc_DTM_50m_avg_ws7_50m FlowLineCurvature_10m
## 2 crosc_ws15_hr_hr Total_Curvature_50m
## 3 slope_DTM_50m_avg_ws3_50m profc_ws15_10m
## 4 Convexity_50m Convergence_Index_hr
## 5 planc_DTM_50m_avg_ws5_50m slope_DTM_50m_avg_ws7_50m
## 6 CrossSectionalCurvature_50m maxic_ws3_hr_hr
## k 5
## 1 minic_ws9_hr_hr
## 2 slope_ws15_10m
## 3 dtm_hr_CONVEX_r30
## 4 profc_DTM_50m_avg_ws11_50m
## 5 minic_ws23_hr_hr
## 6 crosc_ws15_10m
## [1] "OBB error with all predictors of localterrain is 0.472222222222222"
## MDA MDG
## TotalCurvature_10m 0.006871145 1.6716322
## slope_DTM_50m_avg_ws5_50m 0.006275531 1.0674464
## MaximalCurvature_10m 0.004018897 0.9292870
## slope_DTM_50m_avg_ws3_50m 0.003849357 0.9559845
## maxic_ws15_hr_hr 0.003754508 0.8764846
## slope_ws15_10m 0.003743933 1.0338399
## Slope_50m 0.003530568 0.9377814
## slope_ws5_hr_hr 0.003522701 0.8300061
## maxic_DTM_50m_avg_ws5_50m 0.003403529 0.7827860
## maxic_ws7_10m 0.003340688 0.5835607
## [1] "OBB error with all predictors of localterrain is 0.425925925925926"
## [1] "OOB-error: 0.537037037037037 for predictors minic_ws9_hr_hr AND slope_ws19_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 class.error
## 1 11 9 2 2 0.5416667
## 2 10 18 4 2 0.4705882
## 3 3 4 11 8 0.5769231
## 4 0 4 10 10 0.5833333
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.462962962962963"
## [1] "OOB-error: 0.444444444444444 for predictors minic_ws9_hr_hr AND slope_ws19_hr_hr AND minic_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws9_50m AND minic_ws5_hr_hr AND planc_ws9_hr_hr AND profc_DTM_50m_avg_ws3_50m AND crosc_ws13_hr_hr AND dtm_hr_CONVEX_r30 AND planc_ws15_hr_hr AND crosc_DTM_50m_avg_ws11_50m AND crosc_DTM_50m_avg_ws7_50m AND crosc_ws15_hr_hr AND slope_DTM_50m_avg_ws3_50m AND Convexity_50m AND planc_DTM_50m_avg_ws5_50m AND CrossSectionalCurvature_50m AND FlowLineCurvature_10m AND Total_Curvature_50m AND profc_ws15_10m AND Convergence_Index_hr AND slope_DTM_50m_avg_ws7_50m AND maxic_ws3_hr_hr AND slope_ws15_10m AND profc_DTM_50m_avg_ws11_50m AND minic_ws23_hr_hr AND crosc_ws15_10m"
## [1] "confusion OOB"
## 1 2 3 4 class.error
## 1 13 10 0 1 0.4583333
## 2 9 21 3 1 0.3823529
## 3 2 3 15 6 0.4230769
## 4 3 3 7 11 0.5416667
## [1] "1a.4 - Habitat for crops"

## [1] "Prediction error at end is: 0.267099567099567"
## [2] "Prediction error at end is: 0.183982683982684"
## [3] "Prediction error at end is: 0.174458874458874"
## [4] "Prediction error at end is: 0.212121212121212"
## [5] "Prediction error at end is: 0.212121212121212"
## [6] "Prediction error at end is: 0.211688311688312"
## k 1 k 2 k 3
## 1 slope_ws5_hr_hr slope_ws5_10m slope_ws3_hr_hr
## 2 crosc_ws13_hr_hr maxic_ws13_hr_hr longc_ws15_hr_hr
## 3 Slope_hr slope_ws3_10m minic_DTM_50m_avg_ws9_50m
## 4 minic_ws23_hr_hr minic_ws11_10m longc_ws7_10m
## 5 minic_ws5_10m Tangential_Curvature_hr longc_ws3_hr_hr
## 6 planc_ws19_hr_hr maxic_ws7_hr_hr profc_DTM_50m_avg_ws7_50m
## k 4 k 5
## 1 slope_ws3_hr_hr Slope_hr
## 2 minic_ws15_hr_hr longc_ws15_hr_hr
## 3 minic_DTM_50m_avg_ws9_50m slope_ws5_hr_hr
## 4 Convexity_10m profc_ws15_hr_hr
## 5 longc_ws7_hr_hr slope_DTM_50m_avg_ws9_50m
## 6 minic_ws23_hr_hr ProfileCurvature_10m
## [1] "OBB error with all predictors of localterrain is 0.175925925925926"
## MDA MDG
## Slope_hr 0.02615453 2.872642
## slope_ws5_hr_hr 0.02314535 2.777783
## slope_ws3_hr_hr 0.02160712 3.019907
## slope_ws13_hr_hr 0.01695858 2.308600
## slope_ws5_10m 0.01585686 2.070772
## slope_ws11_hr_hr 0.01570075 2.508940
## slope_ws7_10m 0.01493411 1.842763
## slope_ws3_10m 0.01367769 2.098276
## slope_ws29_hr_hr 0.01363032 1.573848
## slope_ws23_hr_hr 0.01261888 1.817747
## [1] "OBB error with all predictors of localterrain is 0.175925925925926"
## [1] "OOB-error: 0.12962962962963 for predictors slope_ws5_hr_hr AND crosc_ws13_hr_hr"
## [1] "confusion OOB"
## 3 4 5 class.error
## 3 5 2 1 0.37500000
## 4 1 58 3 0.06451613
## 5 0 7 31 0.18421053
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.175925925925926"
## [1] "OOB-error: 0.12962962962963 for predictors slope_ws5_hr_hr AND crosc_ws13_hr_hr AND Slope_hr AND minic_ws23_hr_hr AND minic_ws5_10m AND planc_ws19_hr_hr AND slope_ws5_10m AND maxic_ws13_hr_hr AND slope_ws3_10m AND minic_ws11_10m AND Tangential_Curvature_hr AND maxic_ws7_hr_hr AND slope_ws3_hr_hr AND longc_ws15_hr_hr AND minic_DTM_50m_avg_ws9_50m AND longc_ws7_10m AND longc_ws3_hr_hr AND profc_DTM_50m_avg_ws7_50m AND minic_ws15_hr_hr AND Convexity_10m AND longc_ws7_hr_hr AND profc_ws15_hr_hr AND slope_DTM_50m_avg_ws9_50m AND ProfileCurvature_10m"
## [1] "confusion OOB"
## 3 4 5 class.error
## 3 3 5 0 0.6250000
## 4 1 59 2 0.0483871
## 5 0 6 32 0.1578947
## [1] "1c.1 - Average precipitation retention capacity"

## [1] "Prediction error at end is: 0.751082251082251"
## [2] "Prediction error at end is: 0.674458874458874"
## [3] "Prediction error at end is: 0.629437229437229"
## [4] "Prediction error at end is: 0.583982683982684"
## [5] "Prediction error at end is: 0.601298701298701"
## [6] "Prediction error at end is: 0.592207792207792"
## k 1 k 2
## 1 crosc_ws23_hr_hr crosc_DTM_50m_avg_ws9_50m
## 2 planc_ws5_10m profc_ws7_10m
## 3 maxic_DTM_50m_avg_ws3_50m crosc_ws5_10m
## 4 slope_DTM_50m_avg_ws11_50m Convergence_Index_hr
## 5 slope_DTM_50m_avg_ws5_50m planc_DTM_50m_avg_ws9_50m
## 6 dtm_hr_CONVEX_r10 slope_ws15_10m
## k 3 k 4 k 5
## 1 slope_ws29_hr_hr crosc_ws23_hr_hr slope_ws3_hr_hr
## 2 slope_ws5_hr_hr slope_ws11_hr_hr maxic_ws3_10m
## 3 Minimal_Curvature_hr maxic_ws3_10m slope_DTM_50m_avg_ws9_50m
## 4 crosc_DTM_50m_avg_ws11_50m maxic_ws23_hr_hr maxic_ws5_hr_hr
## 5 minic_ws15_10m planc_ws19_hr_hr ProfileCurvature_10m
## 6 TotalCurvature_10m Slope_50m slope_ws15_10m
## [1] "OBB error with all predictors of localterrain is 0.666666666666667"
## MDA MDG
## MaximalCurvature_10m 0.002684007 0.8258478
## maxic_ws3_10m 0.002677553 0.6761451
## TangentialCurvature_10m 0.002452553 0.6587715
## longc_ws3_hr_hr 0.001827930 0.6626584
## slope_DTM_50m_avg_ws11_50m 0.001676560 0.6895404
## slope_DTM_50m_avg_ws5_50m 0.001638520 0.7532804
## Minimal_Curvature_hr 0.001543979 0.5743610
## slope_DTM_50m_avg_ws9_50m 0.001289845 0.8239779
## slope_ws19_hr_hr 0.001244339 0.4637031
## slope_DTM_50m_avg_ws7_50m 0.001162215 0.6094591
## [1] "OBB error with all predictors of localterrain is 0.666666666666667"
## [1] "OOB-error: 0.611111111111111 for predictors crosc_ws23_hr_hr AND planc_ws5_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 5 16 0 3 2 0.8076923
## 2 9 24 3 8 1 0.4666667
## 3 1 6 5 0 0 0.5833333
## 4 2 8 0 8 0 0.5555556
## 5 1 3 1 2 0 1.0000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.648148148148148"
## [1] "OOB-error: 0.574074074074074 for predictors crosc_ws23_hr_hr AND planc_ws5_10m AND maxic_DTM_50m_avg_ws3_50m AND slope_DTM_50m_avg_ws11_50m AND slope_DTM_50m_avg_ws5_50m AND dtm_hr_CONVEX_r10 AND crosc_DTM_50m_avg_ws9_50m AND profc_ws7_10m AND crosc_ws5_10m AND Convergence_Index_hr AND planc_DTM_50m_avg_ws9_50m AND slope_ws15_10m AND slope_ws29_hr_hr AND slope_ws5_hr_hr AND Minimal_Curvature_hr AND crosc_DTM_50m_avg_ws11_50m AND minic_ws15_10m AND TotalCurvature_10m AND slope_ws11_hr_hr AND maxic_ws3_10m AND maxic_ws23_hr_hr AND planc_ws19_hr_hr AND Slope_50m AND slope_ws3_hr_hr AND slope_DTM_50m_avg_ws9_50m AND maxic_ws5_hr_hr AND ProfileCurvature_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 7 14 2 2 1 0.7307692
## 2 8 34 0 3 0 0.2444444
## 3 3 9 0 0 0 1.0000000
## 4 2 11 0 5 0 0.7222222
## 5 3 3 0 1 0 1.0000000
## [1] "1c.1 - Minimum precipitation retention capacity"

## [1] "Prediction error at end is: 0.814718614718615"
## [2] "Prediction error at end is: 0.693073593073593"
## [3] "Prediction error at end is: 0.75021645021645"
## [4] "Prediction error at end is: 0.694372294372294"
## [5] "Prediction error at end is: 0.675324675324675"
## [6] "Prediction error at end is: 0.702597402597403"
## k 1 k 2 k 3
## 1 profc_ws7_10m longc_ws23_hr_hr profc_ws29_hr_hr
## 2 crosc_ws7_hr_hr crosc_ws29_hr_hr maxic_ws15_hr_hr
## 3 planc_ws29_hr_hr minic_ws15_hr_hr minic_ws7_10m
## 4 Total_Curvature_hr crosc_ws5_10m profc_ws3_10m
## 5 crosc_ws3_10m Flow_Line_Curvature_50m minic_ws13_hr_hr
## 6 minic_ws15_hr_hr profc_ws11_10m maxic_ws29_hr_hr
## k 4 k 5
## 1 maxic_ws7_10m planc_ws29_hr_hr
## 2 planc_ws13_hr_hr longc_ws5_10m
## 3 crosc_ws5_10m crosc_ws15_hr_hr
## 4 longc_DTM_50m_avg_ws9_50m crosc_ws15_10m
## 5 CrossSectionalCurvature_hr planc_DTM_50m_avg_ws11_50m
## 6 Tangential_Curvature_50m profc_ws5_10m
## [1] "OBB error with all predictors of localterrain is 0.648148148148148"
## MDA MDG
## TotalCurvature_10m 0.007027602 1.7681303
## Total_Curvature_50m 0.004637173 1.3251695
## maxic_ws15_hr_hr 0.002921044 0.7920725
## minic_ws15_10m 0.002647729 0.8381313
## planc_ws23_hr_hr 0.002493694 0.5586886
## crosc_ws11_10m 0.002244937 0.6849907
## MaximalCurvature_10m 0.002180119 0.6340359
## maxic_ws23_hr_hr 0.002131464 0.5935014
## longc_ws7_hr_hr 0.001609975 0.5133830
## planc_ws29_hr_hr 0.001597839 0.7006263
## [1] "OBB error with all predictors of localterrain is 0.657407407407407"
## [1] "OOB-error: 0.611111111111111 for predictors profc_ws7_10m AND crosc_ws7_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 3 2 3 4 1.0000000
## 2 1 12 5 5 3 0.5384615
## 3 2 3 8 3 6 0.6363636
## 4 2 5 5 11 3 0.5769231
## 5 4 3 2 2 11 0.5000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.694444444444444"
## [1] "OOB-error: 0.5 for predictors profc_ws7_10m AND crosc_ws7_hr_hr AND planc_ws29_hr_hr AND Total_Curvature_hr AND crosc_ws3_10m AND minic_ws15_hr_hr AND longc_ws23_hr_hr AND crosc_ws29_hr_hr AND crosc_ws5_10m AND Flow_Line_Curvature_50m AND profc_ws11_10m AND profc_ws29_hr_hr AND maxic_ws15_hr_hr AND minic_ws7_10m AND profc_ws3_10m AND minic_ws13_hr_hr AND maxic_ws29_hr_hr AND maxic_ws7_10m AND planc_ws13_hr_hr AND longc_DTM_50m_avg_ws9_50m AND CrossSectionalCurvature_hr AND Tangential_Curvature_50m AND longc_ws5_10m AND crosc_ws15_hr_hr AND crosc_ws15_10m AND planc_DTM_50m_avg_ws11_50m AND profc_ws5_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 4 1 5 2 1.0000000
## 2 0 17 2 3 4 0.3461538
## 3 0 2 14 2 4 0.3636364
## 4 1 6 3 11 5 0.5769231
## 5 0 4 0 6 12 0.4545455
## [1] "1c.2 - Retention capacity for heavy precipitation events"

## [1] "Prediction error at end is: 0.471428571428571"
## [2] "Prediction error at end is: 0.352380952380952"
## [3] "Prediction error at end is: 0.287878787878788"
## [4] "Prediction error at end is: 0.296536796536797"
## [5] "Prediction error at end is: 0.296536796536797"
## [6] "Prediction error at end is: 0.296969696969697"
## k 1 k 2
## 1 Minimal_Curvature_50m maxic_ws23_hr_hr
## 2 Longitudinal_Curvature_hr minic_ws29_hr_hr
## 3 longc_DTM_50m_avg_ws3_50m Convexity_50m
## 4 MaximalCurvature_10m slope_DTM_50m_avg_ws9_50m
## 5 maxic_ws13_hr_hr Plan_Curvature_hr
## 6 minic_ws15_10m MaximalCurvature_10m
## k 3 k 4
## 1 crosc_ws7_10m Convexity_50m
## 2 slope_DTM_50m_avg_ws9_50m longc_ws15_10m
## 3 crosc_ws5_hr_hr FlowLineCurvature_10m
## 4 minic_ws5_hr_hr minic_ws3_hr_hr
## 5 profc_DTM_50m_avg_ws5_50m profc_ws11_10m
## 6 Plan_Curvature_50m Flow_Line_Curvature_hr
## k 5
## 1 maxic_ws23_hr_hr
## 2 minic_ws23_hr_hr
## 3 Minimal_Curvature_50m
## 4 maxic_DTM_50m_avg_ws11_50m
## 5 profc_DTM_50m_avg_ws11_50m
## 6 minic_ws15_10m
## [1] "OBB error with all predictors of localterrain is 0.324074074074074"
## MDA MDG
## FlowLineCurvature_10m 0.005909155 1.0104433
## minic_ws7_10m 0.005568158 0.9253621
## crosc_DTM_50m_avg_ws5_50m 0.004630849 0.6995054
## crosc_ws15_10m 0.003865569 0.6760989
## minic_ws29_hr_hr 0.003648695 0.6213939
## Minimal_Curvature_50m 0.003521770 0.6473544
## minic_ws23_hr_hr 0.003410442 0.8414116
## maxic_ws7_10m 0.003320677 0.4325435
## crosc_ws11_10m 0.003207431 0.6300686
## maxic_ws15_hr_hr 0.003193261 0.2923917
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## [1] "OOB-error: 0.222222222222222 for predictors Minimal_Curvature_50m AND Longitudinal_Curvature_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 72 2 1 0 0 0.0400000
## 2 5 5 0 1 0 0.5454545
## 3 6 0 0 0 0 1.0000000
## 4 4 0 0 4 0 0.5000000
## 5 3 2 0 0 3 0.6250000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.314814814814815"
## [1] "OOB-error: 0.25 for predictors Minimal_Curvature_50m AND Longitudinal_Curvature_hr AND longc_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND maxic_ws13_hr_hr AND minic_ws15_10m AND maxic_ws23_hr_hr AND minic_ws29_hr_hr AND Convexity_50m AND slope_DTM_50m_avg_ws9_50m AND Plan_Curvature_hr AND crosc_ws7_10m AND crosc_ws5_hr_hr AND minic_ws5_hr_hr AND profc_DTM_50m_avg_ws5_50m AND Plan_Curvature_50m AND longc_ws15_10m AND FlowLineCurvature_10m AND minic_ws3_hr_hr AND profc_ws11_10m AND Flow_Line_Curvature_hr AND minic_ws23_hr_hr AND maxic_DTM_50m_avg_ws11_50m AND profc_DTM_50m_avg_ws11_50m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 75 0 0 0 0 0.0000000
## 2 7 2 0 2 0 0.8181818
## 3 5 0 0 0 1 1.0000000
## 4 3 1 0 4 0 0.5000000
## 5 6 1 0 1 0 1.0000000
## [1] "1c.3 - groundwater reformation rate"

## [1] "Prediction error at end is: 0.788311688311688"
## [2] "Prediction error at end is: 0.73030303030303"
## [3] "Prediction error at end is: 0.693506493506494"
## [4] "Prediction error at end is: 0.700865800865801"
## [5] "Prediction error at end is: 0.674025974025974"
## [6] "Prediction error at end is: 0.701731601731602"
## k 1 k 2
## 1 slope_ws15_hr_hr LongitudinalCurvature_10m
## 2 profc_ws23_hr_hr maxic_ws7_10m
## 3 maxic_ws19_hr_hr maxic_ws3_hr_hr
## 4 slope_DTM_50m_avg_ws5_50m maxic_ws29_hr_hr
## 5 longc_DTM_50m_avg_ws3_50m slope_DTM_50m_avg_ws5_50m
## 6 minic_ws29_hr_hr GeneralCurvature_10m
## k 3 k 4
## 1 minic_ws15_10m minic_ws11_hr_hr
## 2 FlowLineCurvature_10m maxic_DTM_50m_avg_ws5_50m
## 3 dtm_hr_CONVEX_r5 crosc_ws13_hr_hr
## 4 planc_ws11_10m slope_DTM_50m_avg_ws9_50m
## 5 minic_DTM_50m_avg_ws3_50m planc_ws3_10m
## 6 crosc_ws7_10m maxic_ws9_hr_hr
## k 5
## 1 minic_ws15_10m
## 2 maxic_ws19_hr_hr
## 3 longc_ws3_hr_hr
## 4 maxic_ws15_hr_hr
## 5 maxic_DTM_50m_avg_ws5_50m
## 6 planc_ws7_10m
## [1] "OBB error with all predictors of localterrain is 0.62037037037037"
## MDA MDG
## slope_ws3_10m 0.003745093 1.0041386
## profc_ws7_10m 0.003344335 0.5748368
## slope_ws5_hr_hr 0.002889754 0.7840335
## slope_ws19_hr_hr 0.002721559 0.6605041
## Slope_hr 0.002648683 0.7626642
## slope_ws5_10m 0.002590374 0.6903808
## TotalCurvature_10m 0.002514299 1.3413399
## minic_ws3_hr_hr 0.002448870 0.8295746
## longc_ws3_hr_hr 0.002223869 0.8833354
## slope_ws29_hr_hr 0.002063061 0.6921548
## [1] "OBB error with all predictors of localterrain is 0.657407407407407"
## [1] "OOB-error: 0.601851851851852 for predictors slope_ws15_hr_hr AND profc_ws23_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 1 10 0 2 1 0.9285714
## 2 6 23 2 1 7 0.4102564
## 3 0 5 4 3 2 0.7142857
## 4 1 7 2 6 3 0.6842105
## 5 3 1 3 6 9 0.5909091
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.592592592592593"
## [1] "OOB-error: 0.583333333333333 for predictors slope_ws15_hr_hr AND profc_ws23_hr_hr AND maxic_ws19_hr_hr AND slope_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws3_50m AND minic_ws29_hr_hr AND LongitudinalCurvature_10m AND maxic_ws7_10m AND maxic_ws3_hr_hr AND maxic_ws29_hr_hr AND GeneralCurvature_10m AND minic_ws15_10m AND FlowLineCurvature_10m AND dtm_hr_CONVEX_r5 AND planc_ws11_10m AND minic_DTM_50m_avg_ws3_50m AND crosc_ws7_10m AND minic_ws11_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND crosc_ws13_hr_hr AND slope_DTM_50m_avg_ws9_50m AND planc_ws3_10m AND maxic_ws9_hr_hr AND longc_ws3_hr_hr AND maxic_ws15_hr_hr AND planc_ws7_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 2 9 1 1 1 0.8571429
## 2 1 30 1 1 6 0.2307692
## 3 3 6 1 0 4 0.9285714
## 4 1 6 0 5 7 0.7368421
## 5 2 10 0 3 7 0.6818182
## [1] "1c.4 - Potential for providing nutrients for plants"

## [1] "Prediction error at end is: 0.529004329004329"
## [2] "Prediction error at end is: 0.371428571428571"
## [3] "Prediction error at end is: 0.325108225108225"
## [4] "Prediction error at end is: 0.324675324675325"
## [5] "Prediction error at end is: 0.306493506493506"
## [6] "Prediction error at end is: 0.288311688311688"
## k 1 k 2
## 1 minic_ws5_hr_hr crosc_ws5_hr_hr
## 2 longc_DTM_50m_avg_ws5_50m profc_DTM_50m_avg_ws11_50m
## 3 maxic_ws15_hr_hr planc_DTM_50m_avg_ws11_50m
## 4 profc_ws29_hr_hr minic_ws19_hr_hr
## 5 maxic_ws15_10m Plan_Curvature_50m
## 6 maxic_ws11_10m minic_ws7_hr_hr
## k 3 k 4
## 1 longc_ws29_hr_hr TangentialCurvature_10m
## 2 minic_ws7_10m slope_ws11_10m
## 3 longc_ws7_10m minic_DTM_50m_avg_ws7_50m
## 4 crosc_ws5_hr_hr slope_DTM_50m_avg_ws7_50m
## 5 profc_ws29_hr_hr CrossSectionalCurvature_50m
## 6 CrossSectionalCurvature_50m maxic_DTM_50m_avg_ws5_50m
## k 5
## 1 longc_ws19_hr_hr
## 2 Total_Curvature_50m
## 3 minic_DTM_50m_avg_ws7_50m
## 4 dtm_hr_CONVEX_r30
## 5 planc_DTM_50m_avg_ws11_50m
## 6 minic_ws15_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## MDA MDG
## minic_ws5_hr_hr 0.009522952 1.4110498
## minic_ws7_hr_hr 0.006260161 1.2868091
## minic_ws11_hr_hr 0.006103957 0.9612124
## minic_ws9_hr_hr 0.005741262 1.5157833
## TotalCurvature_10m 0.005000909 1.2009754
## crosc_ws5_hr_hr 0.004862072 1.0021034
## planc_ws5_hr_hr 0.004197141 0.6658763
## minic_ws13_hr_hr 0.004177560 0.9907018
## maxic_ws7_10m 0.003086390 0.6861332
## crosc_ws29_hr_hr 0.003012065 0.4806862
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## [1] "OOB-error: 0.287037037037037 for predictors minic_ws5_hr_hr AND longc_DTM_50m_avg_ws5_50m"
## [1] "confusion OOB"
## 1 3 5 class.error
## 1 54 2 9 0.1692308
## 3 5 3 3 0.7272727
## 5 9 3 20 0.3750000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.296296296296296"
## [1] "OOB-error: 0.25 for predictors minic_ws5_hr_hr AND longc_DTM_50m_avg_ws5_50m AND maxic_ws15_hr_hr AND profc_ws29_hr_hr AND maxic_ws15_10m AND maxic_ws11_10m AND crosc_ws5_hr_hr AND profc_DTM_50m_avg_ws11_50m AND planc_DTM_50m_avg_ws11_50m AND minic_ws19_hr_hr AND Plan_Curvature_50m AND minic_ws7_hr_hr AND longc_ws29_hr_hr AND minic_ws7_10m AND longc_ws7_10m AND CrossSectionalCurvature_50m AND TangentialCurvature_10m AND slope_ws11_10m AND minic_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws7_50m AND maxic_DTM_50m_avg_ws5_50m AND longc_ws19_hr_hr AND Total_Curvature_50m AND dtm_hr_CONVEX_r30 AND minic_ws15_hr_hr"
## [1] "confusion OOB"
## 1 3 5 class.error
## 1 56 2 7 0.1384615
## 3 6 1 4 0.9090909
## 5 8 0 24 0.2500000
## [1] "1c.5 - Potential as a CO2 sink"

## [1] "Prediction error at end is: 0.522943722943723"
## [2] "Prediction error at end is: 0.486580086580087"
## [3] "Prediction error at end is: 0.421212121212121"
## [4] "Prediction error at end is: 0.42987012987013"
## [5] "Prediction error at end is: 0.401731601731602"
## [6] "Prediction error at end is: 0.354112554112554"
## k 1 k 2 k 3
## 1 slope_ws19_hr_hr slope_ws19_hr_hr CrossSectionalCurvature_hr
## 2 PlanCurvature_10m minic_DTM_50m_avg_ws11_50m TotalCurvature_10m
## 3 crosc_ws5_hr_hr planc_ws5_10m FlowLineCurvature_10m
## 4 maxic_ws3_hr_hr planc_ws13_hr_hr planc_DTM_50m_avg_ws7_50m
## 5 TotalCurvature_10m Tangential_Curvature_50m Total_Curvature_50m
## 6 planc_ws3_10m crosc_ws13_hr_hr profc_DTM_50m_avg_ws5_50m
## k 4 k 5
## 1 crosc_ws5_hr_hr crosc_ws19_hr_hr
## 2 slope_ws13_hr_hr profc_ws7_10m
## 3 minic_ws15_hr_hr Tangential_Curvature_50m
## 4 longc_ws7_10m maxic_ws15_10m
## 5 slope_ws15_hr_hr maxic_ws29_hr_hr
## 6 PlanCurvature_10m TangentialCurvature_10m
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## MDA MDG
## TotalCurvature_10m 0.010514204 2.0257496
## maxic_ws15_hr_hr 0.007884619 1.0499995
## MaximalCurvature_10m 0.005602024 1.2579453
## crosc_ws13_hr_hr 0.005087824 0.7986764
## maxic_ws23_hr_hr 0.004971481 0.8665110
## minic_ws5_hr_hr 0.003899458 0.8548583
## slope_ws7_10m 0.003858875 0.7840887
## maxic_ws5_10m 0.003849497 0.6942420
## maxic_ws13_hr_hr 0.003825019 0.9069408
## minic_ws3_hr_hr 0.003466320 0.6909818
## [1] "OBB error with all predictors of localterrain is 0.361111111111111"
## [1] "OOB-error: 0.37962962962963 for predictors slope_ws19_hr_hr AND PlanCurvature_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 43 0 3 5 0 0.1568627
## 2 1 0 0 0 0 1.0000000
## 3 5 0 2 11 0 0.8888889
## 4 4 0 6 22 1 0.3333333
## 5 3 0 0 2 0 1.0000000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.361111111111111"
## [1] "OOB-error: 0.324074074074074 for predictors slope_ws19_hr_hr AND PlanCurvature_10m AND crosc_ws5_hr_hr AND maxic_ws3_hr_hr AND TotalCurvature_10m AND planc_ws3_10m AND minic_DTM_50m_avg_ws11_50m AND planc_ws5_10m AND planc_ws13_hr_hr AND Tangential_Curvature_50m AND crosc_ws13_hr_hr AND CrossSectionalCurvature_hr AND FlowLineCurvature_10m AND planc_DTM_50m_avg_ws7_50m AND Total_Curvature_50m AND profc_DTM_50m_avg_ws5_50m AND slope_ws13_hr_hr AND minic_ws15_hr_hr AND longc_ws7_10m AND slope_ws15_hr_hr AND crosc_ws19_hr_hr AND profc_ws7_10m AND maxic_ws15_10m AND maxic_ws29_hr_hr AND TangentialCurvature_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 44 0 3 4 0 0.1372549
## 2 0 0 1 0 0 1.0000000
## 3 4 0 5 9 0 0.7222222
## 4 6 0 3 24 0 0.2727273
## 5 1 0 0 4 0 1.0000000
## [1] "1d.1 - Potential for retention of heavy metals"

## [1] "Prediction error at end is: 0.638095238095238"
## [2] "Prediction error at end is: 0.5004329004329"
## [3] "Prediction error at end is: 0.435930735930736"
## [4] "Prediction error at end is: 0.427705627705628"
## [5] "Prediction error at end is: 0.408658008658009"
## [6] "Prediction error at end is: 0.399134199134199"
## k 1 k 2
## 1 slope_ws15_10m longc_ws23_hr_hr
## 2 minic_ws7_10m slope_ws29_hr_hr
## 3 dtm_hr_CONVEX_r30 maxic_ws7_10m
## 4 longc_DTM_50m_avg_ws5_50m CrossSectionalCurvature_hr
## 5 planc_DTM_50m_avg_ws5_50m crosc_ws5_hr_hr
## 6 crosc_ws15_hr_hr crosc_ws13_hr_hr
## k 3 k 4
## 1 minic_ws29_hr_hr Total_Curvature_50m
## 2 Flow_Line_Curvature_50m CrossSectionalCurvature_hr
## 3 maxic_ws5_hr_hr minic_ws3_hr_hr
## 4 planc_ws3_hr_hr dtm_hr_CONVEX_r30
## 5 dtm_hr_CONVEX_r20 crosc_DTM_50m_avg_ws7_50m
## 6 FlowLineCurvature_10m TotalCurvature_10m
## k 5
## 1 slope_DTM_50m_avg_ws3_50m
## 2 maxic_ws29_hr_hr
## 3 slope_ws23_hr_hr
## 4 Total_Curvature_hr
## 5 crosc_ws13_hr_hr
## 6 Minimal_Curvature_hr
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## MDA MDG
## minic_ws5_hr_hr 0.010715754 1.6867079
## TotalCurvature_10m 0.006377511 1.2949131
## maxic_ws15_hr_hr 0.004810630 1.0524159
## slope_ws3_10m 0.004442300 0.9790236
## slope_DTM_50m_avg_ws7_50m 0.004344318 0.5440941
## maxic_ws19_hr_hr 0.004108624 0.6347513
## slope_ws5_hr_hr 0.004074841 0.7569981
## slope_ws19_hr_hr 0.003853359 0.7630558
## crosc_ws13_hr_hr 0.003790746 0.8525060
## slope_ws13_hr_hr 0.003670542 0.8060732
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## [1] "OOB-error: 0.37962962962963 for predictors slope_ws15_10m AND minic_ws7_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 34 1 3 0 4 0.1904762
## 2 3 0 0 0 4 1.0000000
## 3 7 0 0 0 3 1.0000000
## 4 3 0 0 0 4 1.0000000
## 5 6 0 2 1 33 0.2142857
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## [1] "OOB-error: 0.324074074074074 for predictors slope_ws15_10m AND minic_ws7_10m AND dtm_hr_CONVEX_r30 AND longc_DTM_50m_avg_ws5_50m AND planc_DTM_50m_avg_ws5_50m AND crosc_ws15_hr_hr AND longc_ws23_hr_hr AND slope_ws29_hr_hr AND maxic_ws7_10m AND CrossSectionalCurvature_hr AND crosc_ws5_hr_hr AND crosc_ws13_hr_hr AND minic_ws29_hr_hr AND Flow_Line_Curvature_50m AND maxic_ws5_hr_hr AND planc_ws3_hr_hr AND dtm_hr_CONVEX_r20 AND FlowLineCurvature_10m AND Total_Curvature_50m AND minic_ws3_hr_hr AND crosc_DTM_50m_avg_ws7_50m AND TotalCurvature_10m AND slope_DTM_50m_avg_ws3_50m AND maxic_ws29_hr_hr AND slope_ws23_hr_hr AND Total_Curvature_hr AND Minimal_Curvature_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 37 0 0 0 5 0.1190476
## 2 4 0 0 0 3 1.0000000
## 3 6 0 2 0 2 0.8000000
## 4 2 0 0 0 5 1.0000000
## 5 8 0 0 0 34 0.1904762
## [1] "1d.2 - Potential for transforming organic contaminants"

## [1] "Prediction error at end is: 0.657142857142857"
## [2] "Prediction error at end is: 0.647186147186147"
## [3] "Prediction error at end is: 0.593073593073593"
## [4] "Prediction error at end is: 0.62034632034632"
## [5] "Prediction error at end is: 0.603463203463203"
## [6] "Prediction error at end is: 0.601298701298701"
## k 1 k 2 k 3
## 1 crosc_ws19_hr_hr MaximalCurvature_10m crosc_ws15_hr_hr
## 2 Minimal_Curvature_hr maxic_ws3_hr_hr crosc_ws19_hr_hr
## 3 maxic_ws29_hr_hr MinimalCurvature_10m maxic_ws15_hr_hr
## 4 minic_DTM_50m_avg_ws9_50m longc_ws15_10m minic_DTM_50m_avg_ws5_50m
## 5 crosc_ws3_hr_hr dtm_hr_CONVEX_r5 longc_DTM_50m_avg_ws9_50m
## 6 slope_DTM_50m_avg_ws3_50m planc_ws15_hr_hr Plan_Curvature_hr
## k 4 k 5
## 1 maxic_ws7_10m crosc_ws13_hr_hr
## 2 maxic_ws23_hr_hr crosc_DTM_50m_avg_ws3_50m
## 3 longc_ws11_10m planc_DTM_50m_avg_ws5_50m
## 4 profc_ws7_hr_hr crosc_ws3_10m
## 5 minic_ws13_hr_hr CrossSectionalCurvature_50m
## 6 longc_ws7_10m minic_ws5_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.62037037037037"
## MDA MDG
## TotalCurvature_10m 0.005572115 1.5331133
## MaximalCurvature_10m 0.004419652 1.2768356
## minic_ws5_hr_hr 0.003956894 1.0044455
## maxic_ws15_hr_hr 0.003554757 1.0823877
## slope_ws7_10m 0.003515566 1.0586123
## slope_ws3_10m 0.002559262 0.6016908
## TangentialCurvature_10m 0.002275492 1.0272114
## slope_ws5_10m 0.002051651 0.5870080
## slope_ws29_hr_hr 0.002026899 0.6302846
## minic_ws7_hr_hr 0.002023120 0.7562407
## [1] "OBB error with all predictors of localterrain is 0.592592592592593"
## [1] "OOB-error: 0.574074074074074 for predictors crosc_ws19_hr_hr AND Minimal_Curvature_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 0 2 4 0 4 1.0000000
## 2 1 8 7 1 5 0.6363636
## 3 3 10 10 0 5 0.6428571
## 4 0 1 1 2 5 0.7777778
## 5 2 6 3 2 26 0.3333333
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.62962962962963"
## [1] "OOB-error: 0.546296296296296 for predictors crosc_ws19_hr_hr AND Minimal_Curvature_hr AND maxic_ws29_hr_hr AND minic_DTM_50m_avg_ws9_50m AND crosc_ws3_hr_hr AND slope_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND maxic_ws3_hr_hr AND MinimalCurvature_10m AND longc_ws15_10m AND dtm_hr_CONVEX_r5 AND planc_ws15_hr_hr AND crosc_ws15_hr_hr AND maxic_ws15_hr_hr AND minic_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws9_50m AND Plan_Curvature_hr AND maxic_ws7_10m AND maxic_ws23_hr_hr AND longc_ws11_10m AND profc_ws7_hr_hr AND minic_ws13_hr_hr AND longc_ws7_10m AND crosc_ws13_hr_hr AND crosc_DTM_50m_avg_ws3_50m AND planc_DTM_50m_avg_ws5_50m AND crosc_ws3_10m AND CrossSectionalCurvature_50m AND minic_ws5_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 2 2 5 0 1 0.8000000
## 2 0 6 10 0 6 0.7272727
## 3 1 9 10 0 8 0.6428571
## 4 0 1 0 0 8 1.0000000
## 5 0 1 7 0 31 0.2051282
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"

## [1] "Prediction error at end is: 0.0554112554112554"
## [2] "Prediction error at end is: 0.0463203463203463"
## [3] "Prediction error at end is: 0.0372294372294372"
## [4] "Prediction error at end is: 0.0372294372294372"
## [5] "Prediction error at end is: 0.0372294372294372"
## [6] "Prediction error at end is: 0.0372294372294372"
## k 1 k 2
## 1 Convexity_10m PlanCurvature_10m
## 2 planc_ws5_10m planc_ws9_hr_hr
## 3 Convexity_50m slope_ws13_hr_hr
## 4 Maximal_Curvature_50m maxic_DTM_50m_avg_ws9_50m
## 5 Minimal_Curvature_50m slope_DTM_50m_avg_ws3_50m
## 6 Plan_Curvature_50m slope_ws15_10m
## k 3 k 4
## 1 CrossSectionalCurvature_50m longc_ws11_10m
## 2 PlanCurvature_10m longc_ws7_10m
## 3 slope_ws13_hr_hr profc_ws5_10m
## 4 Flow_Line_Curvature_50m Minimal_Curvature_50m
## 5 maxic_ws7_hr_hr crosc_ws7_hr_hr
## 6 planc_ws11_hr_hr maxic_ws23_hr_hr
## k 5
## 1 Convexity_10m
## 2 planc_ws9_hr_hr
## 3 PlanCurvature_10m
## 4 minic_ws29_hr_hr
## 5 CrossSectionalCurvature_50m
## 6 minic_ws23_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.037037037037037"
## MDA MDG
## PlanCurvature_10m 0.0013557587 0.39544546
## planc_DTM_50m_avg_ws11_50m 0.0013329326 0.21159295
## longc_DTM_50m_avg_ws5_50m 0.0013044920 0.04848249
## longc_DTM_50m_avg_ws9_50m 0.0011959024 0.07174864
## minic_ws29_hr_hr 0.0010869130 0.05924483
## profc_DTM_50m_avg_ws7_50m 0.0009987909 0.07572986
## planc_ws9_hr_hr 0.0009135926 0.08126357
## longc_ws23_hr_hr 0.0009077445 0.03013816
## longc_ws15_10m 0.0008783968 0.02456702
## profc_DTM_50m_avg_ws9_50m 0.0008627382 0.03370627
## [1] "OBB error with all predictors of localterrain is 0.037037037037037"
## [1] "OOB-error: 0.037037037037037 for predictors Convexity_10m AND planc_ws5_10m"
## [1] "confusion OOB"
## 4 5 class.error
## 4 0 4 1
## 5 0 104 0
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.037037037037037"
## [1] "OOB-error: 0.0462962962962963 for predictors Convexity_10m AND planc_ws5_10m AND Convexity_50m AND Maximal_Curvature_50m AND Minimal_Curvature_50m AND Plan_Curvature_50m AND PlanCurvature_10m AND planc_ws9_hr_hr AND slope_ws13_hr_hr AND maxic_DTM_50m_avg_ws9_50m AND slope_DTM_50m_avg_ws3_50m AND slope_ws15_10m AND CrossSectionalCurvature_50m AND Flow_Line_Curvature_50m AND maxic_ws7_hr_hr AND planc_ws11_hr_hr AND longc_ws11_10m AND longc_ws7_10m AND profc_ws5_10m AND crosc_ws7_hr_hr AND maxic_ws23_hr_hr AND minic_ws29_hr_hr AND minic_ws23_hr_hr"
## [1] "confusion OOB"
## 4 5 class.error
## 4 0 4 1.000000000
## 5 1 103 0.009615385
## [1] "1d.4 - Potential for retention of water-soluble contaminants"

## [1] "Prediction error at end is: 0.795238095238095"
## [2] "Prediction error at end is: 0.600865800865801"
## [3] "Prediction error at end is: 0.572727272727273"
## [4] "Prediction error at end is: 0.554545454545455"
## [5] "Prediction error at end is: 0.526406926406926"
## [6] "Prediction error at end is: 0.535930735930736"
## k 1 k 2
## 1 Minimal_Curvature_hr crosc_ws19_hr_hr
## 2 crosc_ws3_10m CrossSectionalCurvature_hr
## 3 crosc_ws7_hr_hr crosc_ws23_hr_hr
## 4 Slope_50m profc_DTM_50m_avg_ws9_50m
## 5 planc_DTM_50m_avg_ws3_50m crosc_DTM_50m_avg_ws3_50m
## 6 Profile_Curvature_50m profc_ws5_hr_hr
## k 3 k 4
## 1 GeneralCurvature_10m maxic_ws13_hr_hr
## 2 FlowLineCurvature_10m slope_ws29_hr_hr
## 3 crosc_ws3_hr_hr Tangential_Curvature_hr
## 4 crosc_ws15_hr_hr LongitudinalCurvature_10m
## 5 longc_ws3_10m crosc_ws11_hr_hr
## 6 MinimalCurvature_10m longc_DTM_50m_avg_ws3_50m
## k 5
## 1 crosc_ws13_hr_hr
## 2 MinimalCurvature_10m
## 3 crosc_DTM_50m_avg_ws9_50m
## 4 crosc_ws15_hr_hr
## 5 crosc_DTM_50m_avg_ws11_50m
## 6 profc_ws7_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.527777777777778"
## MDA MDG
## slope_ws3_10m 0.003559357 0.8066578
## Slope_hr 0.003275241 0.7043113
## slope_ws23_hr_hr 0.003158708 0.7193676
## minic_ws3_hr_hr 0.003151654 0.9349041
## TangentialCurvature_10m 0.002982387 0.6685046
## CrossSectionalCurvature_hr 0.002937736 0.8851664
## slope_DTM_50m_avg_ws7_50m 0.002867079 0.9348810
## maxic_ws19_hr_hr 0.002682500 0.6667069
## slope_ws11_hr_hr 0.002551860 0.8634671
## slope_ws7_10m 0.002514033 0.5926702
## [1] "OBB error with all predictors of localterrain is 0.5"
## [1] "OOB-error: 0.5 for predictors Minimal_Curvature_hr AND crosc_ws3_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 33 2 4 6 2 0.2978723
## 2 9 0 3 2 0 1.0000000
## 3 7 0 5 4 1 0.7058824
## 4 6 1 1 14 0 0.3636364
## 5 3 1 2 0 2 0.7500000
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.490740740740741"
## [1] "OOB-error: 0.5 for predictors Minimal_Curvature_hr AND crosc_ws3_10m AND crosc_ws7_hr_hr AND Slope_50m AND planc_DTM_50m_avg_ws3_50m AND Profile_Curvature_50m AND crosc_ws19_hr_hr AND CrossSectionalCurvature_hr AND crosc_ws23_hr_hr AND profc_DTM_50m_avg_ws9_50m AND crosc_DTM_50m_avg_ws3_50m AND profc_ws5_hr_hr AND GeneralCurvature_10m AND FlowLineCurvature_10m AND crosc_ws3_hr_hr AND crosc_ws15_hr_hr AND longc_ws3_10m AND MinimalCurvature_10m AND maxic_ws13_hr_hr AND slope_ws29_hr_hr AND Tangential_Curvature_hr AND LongitudinalCurvature_10m AND crosc_ws11_hr_hr AND longc_DTM_50m_avg_ws3_50m AND crosc_ws13_hr_hr AND crosc_DTM_50m_avg_ws9_50m AND crosc_DTM_50m_avg_ws11_50m AND profc_ws7_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 40 0 3 4 0 0.1489362
## 2 9 0 2 3 0 1.0000000
## 3 7 1 2 7 0 0.8823529
## 4 6 0 4 12 0 0.4545455
## 5 1 0 2 5 0 1.0000000
## [1] "1d.5 - Potential as buffer for acidic contaminants"

## [1] "Prediction error at end is: 0.759307359307359"
## [2] "Prediction error at end is: 0.684848484848485"
## [3] "Prediction error at end is: 0.674458874458874"
## [4] "Prediction error at end is: 0.674891774891775"
## [5] "Prediction error at end is: 0.637662337662338"
## [6] "Prediction error at end is: 0.647186147186147"
## k 1 k 2 k 3
## 1 crosc_ws29_hr_hr crosc_ws23_hr_hr planc_ws7_hr_hr
## 2 FlowLineCurvature_10m planc_ws11_10m maxic_ws23_hr_hr
## 3 planc_DTM_50m_avg_ws11_50m longc_ws15_10m planc_DTM_50m_avg_ws3_50m
## 4 crosc_ws3_hr_hr planc_ws3_10m crosc_ws19_hr_hr
## 5 crosc_ws13_hr_hr profc_ws15_10m crosc_ws5_hr_hr
## 6 planc_ws11_10m crosc_ws15_10m planc_ws23_hr_hr
## k 4 k 5
## 1 minic_DTM_50m_avg_ws9_50m slope_ws3_10m
## 2 crosc_ws5_hr_hr planc_ws5_hr_hr
## 3 slope_DTM_50m_avg_ws9_50m crosc_ws19_hr_hr
## 4 Slope_50m crosc_ws5_hr_hr
## 5 longc_ws5_hr_hr maxic_ws11_10m
## 6 longc_ws15_hr_hr slope_ws13_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.583333333333333"
## MDA MDG
## maxic_DTM_50m_avg_ws11_50m 0.004756981 1.2730419
## crosc_ws5_hr_hr 0.003770875 0.7768003
## minic_ws11_hr_hr 0.003744165 0.8244634
## maxic_DTM_50m_avg_ws9_50m 0.003010718 0.7370658
## slope_ws13_hr_hr 0.002727887 0.6191520
## longc_DTM_50m_avg_ws11_50m 0.002544911 0.7190560
## planc_ws23_hr_hr 0.002420265 0.7082539
## Total_Curvature_50m 0.002318242 1.0620415
## profc_DTM_50m_avg_ws11_50m 0.002313833 0.6941452
## slope_ws23_hr_hr 0.001891837 0.5294058
## [1] "OBB error with all predictors of localterrain is 0.574074074074074"
## [1] "OOB-error: 0.62962962962963 for predictors crosc_ws29_hr_hr AND FlowLineCurvature_10m"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 25 6 7 0 1 0.3589744
## 2 10 6 4 3 1 0.7500000
## 3 7 2 3 6 2 0.8500000
## 4 2 3 8 3 0 0.8125000
## 5 3 0 2 1 3 0.6666667
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.574074074074074"
## [1] "OOB-error: 0.546296296296296 for predictors crosc_ws29_hr_hr AND FlowLineCurvature_10m AND planc_DTM_50m_avg_ws11_50m AND crosc_ws3_hr_hr AND crosc_ws13_hr_hr AND planc_ws11_10m AND crosc_ws23_hr_hr AND longc_ws15_10m AND planc_ws3_10m AND profc_ws15_10m AND crosc_ws15_10m AND planc_ws7_hr_hr AND maxic_ws23_hr_hr AND planc_DTM_50m_avg_ws3_50m AND crosc_ws19_hr_hr AND crosc_ws5_hr_hr AND planc_ws23_hr_hr AND minic_DTM_50m_avg_ws9_50m AND slope_DTM_50m_avg_ws9_50m AND Slope_50m AND longc_ws5_hr_hr AND longc_ws15_hr_hr AND slope_ws3_10m AND planc_ws5_hr_hr AND maxic_ws11_10m AND slope_ws13_hr_hr"
## [1] "confusion OOB"
## 1 2 3 4 5 class.error
## 1 29 6 2 1 1 0.2564103
## 2 10 8 3 2 1 0.6666667
## 3 6 4 6 3 1 0.7000000
## 4 4 2 2 6 2 0.6250000
## 5 2 1 3 3 0 1.0000000